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A  branch-and-bound  solution  algorithm  and  a  computer  program 
implementing  this  algorithm  are  developed  to  solve  a  multiactivity 
multifacility  capacity  constrained  0-1  assignment  problem.  The  math¬ 
ematical  formulation  for  such  a  problem,  called  problem  (P),  is  to  find 

x..  and  y,  values  that: 
ij  -'k 

m  n  p 

Minimize  E  E  a . .  x .  .  +  E  b.  y.  (i) 

i-1  J-l  «  1J  k-1  k  k 
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where  i  ,  j  ,  k  are  indices  for  designs,  activities,  and  facilities, 
respectively;  x^  has  value  1  if  and  only  if  activity  j  uses  design 

i  ,  and  y^  has  value  1  if  and  only  if  facility  k  is  used.  A  design 

involves  the  use  of  one  or  more  facilities,  and  the  same  design  may  be  used 
by  several  activities. 


Problem  (P)  has  the  objective  of  minimizing  the  sum  of  a..'s  — 

the  variable  costs  due  to  the  assignments  of  activities  to  designs,  and 
b^'s  —  the  fixed  costs  due  to  the  facilities  used.  Constraints  (ii) 


and  (iv)  ensure  that  each  activity  is  assigned  to  c  single  design. 

d.  is  the  capacity  required  at  facility  k  if  activity  j  uses 
1  j  k 


Each 

design 


i  ,  and  is  thus  equal  to  zero  if  design  i  does  not  involve  the  use  of 
facility  k  .  Constraints  (iii),  therefore,  ensure  that  for  each  facility 
k  used,  the  total  capacity  required  does  not  exceed  the  capacity  available 
s^  .  The  difficulty  in  solving  problem  (P)  stems  from  the  indirect  relation¬ 


ship  between  the  assignments  and  facilities,  i.e.,  an  assignment 
bears  on  all  the  constraints  (iii)  for  which  d  ^  is  positive, 
therefore,  on  several  y^  variables. 


Aij 

and , 
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The  branch-and-bound  solution  algorithm  uses  Lagrangian  relaxation 
as  a  basic  step  in  obtaining  lower  bounds.  In  addition,  it  includes  several 
operational  rules,  such  as  a  branching  rule  for  a  judicious  choice  of  the 
branching  variable,  a  capacity  rule  to  eliminate  infeasible  assignments, 
and  a  bounding  rule  to  eliminate  non-optimal  assignments. 


This  dissertation  includes  relevant  background  leading  to  the 
formulation  of  problem  (P) ,  mathematical  development  of  the  branch-and- 
bound  solution  algorithm,  a  detailed  test  example,  and  computational  test 
results  using  the  computer  program.  The  areas  of  application  are  identified, 
and  suggestions  for  further  improvement  of  the  branch-and-bound  solution 
algorithm  are  included. 

The  computer  program  has  been  written  in  FORTRAN  IV.  A  detailed 
description  of  the  computer  program  and  guidelines  for  its  use  are  included 
in  a  separate  document  entitled  "Program  Description  and  User's  Guide  for 
ZIPCAP — a  Zero-one  Integer  Program  to  solve  multiactivity  multifacility 
Capacity-constrained  Assignment  Problems."  Although  developed  for  capac¬ 
itated  problems,  the  computer  program  can  also  be  used  to  solve  uncapacitated 
problems  in  which  it  is  assumed  that  the  facilities  have  infinite  capacity. 
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1.  INTRODUCTION 


Multiact ivity  multifacility  assignment  problems  arise  in  such 
diverse  areas  as  public  health  care  systems  and  private  multi-echelon 
inventory/distribution  systems.  Such  systems  involve  the  assignment  of 
activities  or  tasks  to  groups  of  facilities  in  such  a  way  that  total 
system  cost  is  minimized.  The  total  system  cost  has  components  (fixed 
costs)  that  depend  on  the  facilities  actually  used  as  well  as  components 
(variable  costs)  that  depend  solely  on  the  assignment  made.  Most  recently 
[Gross,  Pinkus,  and  Soland  (1979)]  tiiere  has  been  interest  in  including 
facility  capacity  constraints  as  well.  For  this  kind  of  problem,  i.e., 
a  multiactivity  multifacility  capacity-constrained  0-1  assignment  problem, 
we  have  developed  a  solution  algorithm  of  the  branch-and-bound  type  and 
a  computer  program  based  upon  it. 

The  computer  program  and  guidelines  for  its  use  are  described  in 
a  separate  document  (Chhabra  and  Soland  (1980)]  titled  "Program  Descrip¬ 
tion  and  User's  Guide  for  ZIPCAP  —  a  Zero-one  Integer  Program  to  solve 
multiactivity  multifacility  Capacity-constrained  Assignment  Problems." 

This  document  describes  the  development  of  the  solution  algorithm 
and  computational  test  results  using  the  computer  program.  Suggestions 
for  further  improvement  in  the  solution  algorithm  are  also  included. 

This  chapter  reviews  the  relevant  literature,  provides  background 
leading  to  the  mathematical  formulation  of  the  multiactivity  mul t ifaci 1 i ty 
capacity-constrained  0-1  assignment  problem,  called  problem  (P) ,  and 
includes  potential  areas  of  application.  The  theoretical  base  for  develop¬ 
ing  the  algorithm/methodology  are  described  in  Chapter  2.  Various  components 
of  the  methodology  are  covered  in  detail  in  Chapter  3.  Chapter  4  provides 
an  overview  of  the  computational  procedure  and  the  computer  program,  whereas 
computational  test  results  are  given  in  Chapter  5.  Suggestions  for  further 
research  and  potential  improvements  in  the  algorithm  are  included  in  Chapter 
6. 
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It  may  be  noted  that  the  basic  terminology,  described  below, 
in  the  formulation  of  problem  (P)  includes:  activities  that  must  be 
assigned,  facilities  which  serve  the  activities,  designs  involving  one 
or  more  facilities,  fixed  costs  associated  with  tne  facilities,  and 
variable  costs  associated  with  the  assignment  of  activities  to  designs. 

The  following  review  of  the  relevant  literature  starts  with  the 
classical  assignment  problem  and  leads  to  the  fonnulation  of  problem  (P). 
Different  authors  have  used  various  terminologies  in  describing  relevant 
formulations.  In  the  following  discussion,  the  original  terminologies 
are  used,  and  are  followed  by  our  equivalent  terminology,  where  appro¬ 
priate,  shown  in  parenthesis. 

1 . 1  Generalized  Assignment  Problem 

In  a  classical  assignment  problem  [Hillier  and  Lieberman  (1980)], 
the  purpose  is  to  find  optimal  pairs  of  agents  ard  tasks  or  activities. 
Each  task  is  assigned  to  a  single  agent,  and  each  agent  is  given  a 
single  task,  and  the  suitability  of  a  particular  set  of  assignments 
is  determined  by  a  single  criterion  function  such  as  minimization  of 
cost.  In  a  generalized  assignment  problem  (GAP),  several  tasks  can  be 
assigned  each  agent,  subject  to  the  resources  available  to  the  various 
agents  [Ross  and  Soland  (1975)],  e.g.,  assigning  software  development 
tasks  to  programmers  and  assigning  jobs  to  computers  in  a  computer 
network. 

A  variety  of  well-known  facility  location  and  location-alloca¬ 
tion  problems  have  been  shown  to  be  equivalent  to,  and  therefore 
solvable  as  GAP's  [Ross  and  Soland  (1977)].  Here,  in  general,  the 
tasks  represent  demand  centers  for  a  good  or  service,  and  the  agents 
represent  supply  centers  to  be  established  at  potential  sites  or 
locations.  Each  demand  center  must  be  supplied  from  a  supply  center. 

A  fixed  cost  is  incurred  for  each  supply  center  established,  and,  in 
addition,  there  is  a  cost  incurred  for  each  unit  processed  at  a  supply 
center  and  transportation  costs  incurred  for  the  units  sent  from  supply 
centers  to  demand  centers.  The  problem  may  be  "uncopacitated"  —  when 
there  is  no  limit  to  the  number  of  units  that  may  be  processed  by 


a  supply  a'HUT,  or  "capacitated"  —  when  there  are  restrictions  on  the 
number  of  units  that  may  be  processed.  The  objective  is  to  select  supply 
center  locations  and  set  up  a  distribution  assignment  so  that  the  total 
cost  is  minimized. 

1 . 2  Mult iact ivity  Mul t iiac il ity  Uncapacitated 
Assignment  Problem 

A  salient  feature  of  the  above  facility  location  problems  is  that 
each  demand  tenter  (activity)  is  assigned  to  a  single  supply  center 
(facility).  Sometimes,  however,  it  may  be  desirable  to  assign  an  activity 
to  more  than  one  facility.  This  leads  to  the  concept  of  "design,"  and 
the  multiactivity  multifacility  assignment  problem  [Pinkus,  Gross,  and 
Soland  (197J) J .  Before  describing  such  a  problem,  some  terminology  is 
considered  first. 

A  design  involves  the  use  of  one  or  more  facilities,  and  represents 
a  meaningful  configuration  of  facilities  along  with  a  meaningful  strategy 
for  using  them  —  as  illustrated  in  the  following  examples. 

Consider  five  facilities  and  their  locations  as  shown  in  Figure 
1(a).  (From  practical  considerations,  these  may  be  existing  and/or 
potential  locations.)  Three  of  the  possible  designs  are  shown  in 
Figures  1(b)  to  1(d).  Design  1  is  completely  cent-alized  since  it  uses 
only  one  facility,  whereas  design  3  is  completely  decentralized  since  it 
uses  ail  the  facilities. 
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(a)  Design  2:  three  facili¬ 
ties — locations  0.  0 
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Figure  1. 
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(d)  Design  3:  five  facilities 
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Examples  of  alternate  designs 
for  a  system  of  five  facilities 
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It  is  possible  for  several  designs  to  have  the  same  facilities 
but  different  configuration  and  strategies  for  using  these  facilities, 
e.g.,  a  multiproduct  multi-echelon  inventory  system  [Gross,  Pinkus, 
and  Soland  (1979)].  Figure  2(a)  shows  design  1  containing  certain 
facilities  (warehouses)  at  the  central,  regional,  and  local  levels 
or  echelons.  Figure  2(b)  shows  design  2  with  the  same  facilities 
but  having  a  different  configuration. 
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(a)  Design  1  (h)  Design  2 
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Figure  2. 


Example  of  alternative  designs  having 
the  sane  facilities  but  different 
configuration 


The  distribution  of  a  given  activity  at  various  facilities  under  design  1 
would  be  different  than  under  design  2,  depending,  of  course,  on  the 
inventory  policies.  This  results  in  different  variable  costs  (described 
later)  for  that  activity  under  design  1  as  against  design  2.  In  fact, 
it  is  possible  to  have  a  situation  where  two  or  more  designs  have  the 
same  facilities  and  the  same  configuration  but  different  strategies, 
resulting  in  different  variable  costs.  For  example,  one  strategy  might 
specify  an  equal  distribution  of  a  specific  activity  over  the  various 
facilities,  whereas  another  strategy  could  impose  a  different  distribu¬ 
tion  scheme  over  the  same  facilities. 

In  general,  if  a  system  is  to  be  composed  of  at  most  p 

facilities,  the  number  of  alternative  designs  is  2^-1  if  no  two  designs 
have  the  same  facilities.  However,  with  the  same  facilities  but  different 
configurations  and  strategies,  the  number  of  alternative  designs  could  be 
much  higher.  In  practice,  it  is  possible  to  eliminate  a  majority  of 
alternative  designs  because  of  geographical,  political,  economical,  and 
other  factors. 
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i'he  mnl t iactiv ity  mu  1 1  if  ac  i  L  ity  assignment  problem  seeks 
minimi/,  at ion  of  some  measure  of  total  system  cost  such  as,  total 
expected  cost  over  a  given  time  period  or  total  discounted  cost  over 
the  lifetime  of  t lie  system.  The  system  cost  will  include  investment 
costs  for  building  or  leasing  the  system,  operating  costs  for  operation 
and  maintenance  of  the  system,  and  the  costs  for  providing  necessary 
services.  Both  the  investment  costs  and  the  opeiating  costs  have 
fixed  as  well  as  variable  components  [Ross  and  Soland  (1980)].  The 
fixed  components  include  those  costs  associated  with  the  facilities 
ot  j  given  design  which  are  independent  of  the  activities  served. 

Such  costs  are  called  fixed  costs.  On  the  other  hand,  the  variable 
components  and  the  service  costs  include  those  costs  which  are 
completely  dependent  on  the  service  demand  of  the  activities  at 
the  various  facilities  in  a  given  design.  Such  costs  are  called 
variable  costs.  By  definition,  both  the  fixed  costs  and  the  variable 
costs  are  relative  terms. 

An  equivalent  formulation  of  the  multiact ivity  multifacility 
assignment  problem  defined  by  Pinkus,  Gross,  ant.  Soland  (1973)  is  as 
follows. 


Let  a  .  . 
ij 


b 


k 


b 


ik 


variable  cost  of  activity  j  using  design  i 
(i-1, . . . ,m;  j— l,..,,n) 

fixed  cost  of  facility  k  (k=l,...,p) 

1  if  facility  k  is  included  in  design  i  , 


0  otherwise. 


is  defined  as: 
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The  decision  variable  x. . 

ij 


x _  =  1  if  activity  j  uses  design  i  , 
=  0  otherwise. 


Then,  the  uncapacitated  assignment  problem  called  problem 
(PU)  is  to  find  x  values  that: 


(PU)-C 


Minimize 


subject  to 


where 


m  n 

Z  I  a.,  x..  +  Z  b,  ul  £  b.,  Z  x 

1J  1J  k=l  k  li=l  lk  j-1  « 


i=l  J-1 


m 


Z  x  .  .  =  1 


i=l 
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for  j-1, ... ,n 


x..  =  0  or  1  for  all  i  and  j 
ij  J 


u(-)  =  0  if  (•)  <  0  , 
=  1  if  (•)  >  0  . 


(1) 

(2) 

(3) 


The  objective  function  of  this  problem  coisists  of  two  distinct 
parts.  The  first  part  represents  the  total  variable  cost,  and  the 
second,  the  total  fixed  cost  of  the  system.  Constraints  (2)  and  (3) 
ensure  that  each  activity  is  assigned  to  a  single  design.  Of  course, 
the  optimal  solution  may  involve  the  use  of  more  chan  one  design. 

Problem  (PU)  is  a  0-1  nonlinear  programming  problem  (because  of 
the  step  function  u  ),  and  a  branch-and-bound  algorithm  using  linear 
underestimates  for  the  nonlinear  part  of  the  objective  function  has  been 
described  in  Pinkus,  Gross,  and  Soland  (1973).  A  heuristic  procedure 
for  this  problem  is  given  by  Khumawala  and  Stinson  (1980)  in  an  unpublished 
paper.  This  procedure  is  an  extension  of  some  earlier  work  [Khumawala 
(1973) ] . 

1 . 3  Adding  Cap  icity  Constraints  —  P ro hi  cm  (P) 

A  weakness  of  problem  (PU)  is  that  it  assumes  unlimited  capacity 
available  at  each  facility  in  terms  of  the  activities  using  a  given 
facility.  In  practice,  a  facility  may  not  have  the  capability  to 
serve  every  activity,  and  may  have  restrictions  as  to  the  total 
capacity  available  to  handle  more  than  one  activity. 
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Let  Sj  ~  capac i t v  ova  i  1  ahl  e  at  faciliiv  U  ,  and 

d..^  =  capacity  required  at  facility  k  for  activity  j 
when  activity  j  uses  design  i  . 

If  design  i  does  not  include  facility  k  ,  then  d  =  0  for  all  j 

1 1  k 

Define  the  decision  variable  y,  as: 

J  k 

y^  =  1  if  facility  k  is  used, 

=  0  otherwise  . 

Then  the  assignment  problem  (Gross,  Pinkus,  and  Soland  (1979)], 

called  problem  (P)  is  to  find  x..  and  y,  values  that: 

✓  ij  ■’k 


Minimize  I  E  a.,  x.  +  E  b,  y, 

i=l  j=i  ^  k=1  k  k 


subject  to 


E  x. .  =  1 


j=l . n 


if1  dijk  xij  1  Vk  k=1--*’p 

xij’  yk  =  0  or  1  for  a11  i>j»k 


Constraints  (5)  of  problem  (F)  ensure  that  the  capacities  available  at 
the  facilities  are  not  violated.  Problem  (P)  is,  thus,  a  multiactivity 
multifacility  capacity-constrained  0-1  assignment  problem,  as  compared 
to  problem  (PU)  which  is  uncapacitated.  In  problem  (P) ,  constraints  (2) 
along  with  the  part  of  constraints  (6)  involving  the  x^'s  ensure  that 

each  activity  is  assigned  to  a  single  design.  Of  course,  the  optimal 
solution  may  result  in  the  use  of  more  than  one  design. 

For  an  example  of  five,  facilities  and  three  designs  as  shown 

in  Figures  1(b)  to  1(d),  and  four  activities;  the  matrix  [a  .  lb,  Id...] 

ij  1  k1  ljk 

is  as  shown  in  Figure  3. 

1.3. L  Comparison  with  the  uncapacitated  assignment  problem. 

Comparison  of  the  capacitated  problem  (P)  with  the  uncapacitated 
problem  (PU)  shows  that  the  objective  functions  (1)  and  (4)  are  equiva¬ 
lent  and  constraints  (2)  in  each  are  the  same.  Constraints  (5)  serve 
to  impose  the  capacity  constraints  and  at  the  same  time,  for  a  given 
design,  the  relevant  facilities  are  forced  in  the  solution.  For  an 


,  -.i  ^ - 
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x  equal  to  1  ,  all  the  facilities  with  d..,  >  0  .nust  have  y, 
ij  n  ljk  7k 

values  equal  to  1  in  order  to  satisfy  (5)  and  the  corresponding  fixed 

costs  are  therefore  included  in  (4).  If  y^_  =  0  and  d„^  >  ® 

then  x..  must  be  0  in  order  to  satisfy  (5). 


ProbLem  (P)  lias  been  formulated  as  a  0-1  linear  programming 
problem  whereas  problem  (PLi)  was  formulated  as  a  0-1  nonlinear 
programming  problem. 


case  of 
design 


Note  that  problem  (PU)  can  be  easily  obtained  as  a  special 

problem  (P)  by  letting  d..,  equal  ]  (for  all  j  )  if 

ij  k 

i  uses  facility  k  ,  and  setting  all  s^  equal  to  n  . 


In 


other  words,  the  corresponding  formulation  is  to  find 
values  that: 


x .  . 
ij 


and 


y 


k 


(PI) 


m  n 


Minimize  L  L  a..x..  +  Z  b,  y, 

■  i  •  .  ij  ij  .  .  k  k 

i=l  j=l  J  J  k=l 


m 

l 

i=l 


subject  to  £  x  =  1 


j=l, . . . ,n 


m  n 

Z  e  ..  Z  x  .  .  <  ny.  k=l ,  .  .  .  .p 

lk  .  ,  lj  —  k 


1=1  j=l 


x..,y,  =  0  or  1  for  all  i,j,k 
ij’  k 


I  where 


'ik 


1  if  design  i  uses  facility  k  , 
0  otherwise  . 


(4) 


(2) 

(7) 

(6) 


1.3.2  Comparison  with  the  fixed-charge  location-allocation  problem. 

Problem  (P)  bears  a  resemblance  to  the  well-known  fixed-charge 
location-allocation  problem  or  capacitated  facility  location  problem 
[Geoff rion  (1975);  Ross  and  Soland  (1977)].  There  are,  however, 
very  significant  differences  between  the  two.  In  order  to  point  out 
these  differences,  here  is  a  statement  of  the  location-allocation 
problem  (Ij\)  as  given  by  Gross,  Pinkus,  and  Soland  (1979)  in  a  form 


_  10 


similar  to  that  oil  problem  (P) . 


Find  x  and  y^  values  that 


(LA) 


Minimize 


subject  to 


p  n  p 

A  A  x“j  +  A ^ 


£  X  =  1 

k-1  kJ 


A  di  xkj  -  v* 


j~i»  * • • 


k=l , . , 


Xkj  -  0  ’  yk  '  0  °r  1  '-or  a11  j 

and  k 


(8) 


(9) 


(10) 


(11) 


Here  x^  represents  the  fraction  of  customer  (activity)  j's  demand 
that  is  supplied  by  a  facility  at  location  k. 

The  most  important  distinction  between  problen  (LA)  and  problem 
(P)  is  the  relationship  between  assignments  and  facilities.  In  problem 
(LA)  there  is  a  direct  connection  between  the  assignments  made  and  the 
facilities  required,  and  each  assignment  affects  cnly  one  facility,  i.e., 
the  assignment  >  0  has  a  bearing  on  only  one  of  the  constraints 


(10)  and,  therefore,  on  only  one  variable  y 


On  the  other  hand , 


in  problem  (P),  the  connection  between  the  assignments  made  and  the 
facilities  required  is  indirect,  and  each  assignment  can  affect  several 
facilities,  i  e.,  the  assignment  x„  =  1  bears  on  all  of  the  constraints 

(5)  for  which  d.  >  0  and,  therefore,  on  several  variables  y,  . 

1 J  K.  K, 

Another  distinction  is  the  relative  difficulty  of  the  two  problems. 
While  problem  (LA)  is  not  easy  to  solve,  branch-and-bound  approaches 
have  been  successful  in  dealing  with  it  because  oncz  values  are  specified 
for  the  y^  ,  the  x.^  are  found  by  solving  a  transportation  problem. 

Problem  (LA)  becomes  more  difficult  if  the  constraints  x  >  0  in  (10) 

kj  - 

are  replaced  by  x  =  0  or  1  in  order  to  preclude  supply  of  customer 
(activity)  j's  demand  by  more  than  one  facility.  With  this  change. 
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problem  (LA)  may  be  treated  as  a  generalized  assignment  problem  and 
is  solvable  using  an  efficient  branch-and-bound  algorithm  [Ross  and 
Soland  (1977)].  Problem  (P)  is  more  difficult  than  this  variation 
of  problem  (LA)  because  of  the  above  stated  indirect  connection 
between  the  assignments  and  the  facilities.  Even  after  values  have 
been  specified  for  all  the  y^  ,  problem  (P)  remains  a  difficult 

0-1  linear  programming  problem  because  of  the  interaction  of  the 
constraints. 

1.3.3  Solving  problem  (P). 

The  capacitated  problem  (P)  has  mn+p  0-1  variables  and 

n+p  constraints,  so  the  problem  dimensions  may  be  large  from  practical 
considerations.  For  example,  with  m=n=30  and  p=20  ,  problem  (P) 
has  920  variables  and  50  constraints.  The  0-1  LP  computer  codes 
generally  available  are  limited  in  terms  of  problem  size.  For  example, 
the  code  used  by  Gross,  Pinkus,  and  Soland  (1979)  can  handle  up  to  40 
variables  and  20  constraints.  A  better  and  more  efficient  code 
[Geoffrion  and  Nelson  (1968)]  allows  up  to  90  variables  and  50 
constraints.  This  fact,  together  with  the  structure  of  problem  (P) 
suggests  that  a  specialized  algorithm  could  be  developed  that  would 
be  more  efficient  for  practical  problems  than  the  general  integer 
linear  programming  algorithms  (on  which  the  available  codes  are 
based)  . 

With  the  above  background  in  mind,  the  development  of  the 
solution  algoritlim  and  the  computer  program  to  solve  problem  (P) 
was  undertaken  and  is  described  in  Chapters  2  thrsugh  4. 

1 . 4  Areas  of  App 1 ication 

The  solution  algorithm  and  the  computer  program  are  designed 
to  solve  a  multiactivity  multifacility  capacity-constrained  0-1 
assignment  problem,  i.e.,  one  which  can  be  formulated  as  problem  (P) . 

The  basic  elements  of  such  a  problem  are  activities  that  must 
be  assigned,  facilities  and  their  meaningful  configurations 


represented  as  designs,  the  fixed  and  variable  costs,  and  the 
capacity  requirements  of  the  activities. 


The  formulation  (P)  applies  to  existing  and/or  proposed 
facilities.  In  other  words,  it  is  useful  for  a  situation  where  the 
decision  may  be  to  delete  some  of  the  existing  facilities,  as  well  as 
for  a  situation  where  the  decision  involves  a  selection  out  of  a  set 
of  proposed  facilities. 

Table  1  includes  examples  of  areas  where  formulation  (P)  is 
applicable.  Within  each  application  area,  activities  and  facilities 
are  identified.  The  implications  of  designs,  variable  costs,  and 
fixed  costs  are  apparent. 

Obtaining  the  values  of  the  data  elements  b, ,  d  , ,  s  , 

k  ij  k  k 

and  in  particular  a  ,  can  be  a  simple  or  a  complex  exercise  depending 

on  tl'.e  particulars  of  the  application,  and  the  nature  of  the  components 
comprising  these  elements.  For  example,  in  designing  multi-echelon 
inventory  systems  [Gross,  Pinkus,  and  Soland  (1979)],  a_  represents 

the  inventory  cost  of  product  (activity)  j  using  echelon  structure 
(design)  i  and  b^  represents  the  fixed  cost  of  installation  (facility) 

k  .  The  inventory  cost  a_  includes  the  cost  of  procurement,  carrying 

inventory,  filling  orders,  and  stockouts.  The  value  a^  ,  and  associated 

inventory  stockage  policies,  are  arrived  at  by  solving  a  multi-echelon 
inventory  problem.  In  other  words,  for  product  j  stocked  under  echelon 


structure  i 
installation 


,  optimal  inventory  policies  are  determined,  at  each 
of  the  structure,  which  yields  a^  .  The  facility  fixed 


cost  b^  includes  the  capital  expenditure  for  building  the  installation. 


along  with  a  number  of  fixed  costs  associated  with  operating  it,  such  as 
administrative  expenses,  the  expense  of  renting  the  facility  (if  it  is 
not  built),  .nd  certain  other  fixed  operating  expenses. 


In  the  case  of  designing  a  support  system  for  repairable 


EXAMPLES  op  application  apeas 


*Gross  and  Pinkus  (1979) 
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items  [Gross  and  Pinkus  (1979)],  a^  represents  the  total  variable 

cost  if  unit  type  (activity)  j  is  repaired  under  design  i  .  The 
set  of  parameters  taken  into  consideration  to  compute  this  cost  for 
each  unit  type  includes  such  things  as  varying  population  sizes, 
failure  rates,  average  repair  times,  costs  associated  with  their 
repair,  the  purchase  and  storage  of  spares,  the  purchase  of  repair 
channels,  and  travel  to  depots  (facilities)  for  repair.  A  computer 
program  is  used  to  solve  a  spares  and  server  provisioning  problem,  and 
the  results  provide  the  basic  information  to  compute  a^  . 

Thus,  in  general,  the  data  elements  of  problem  (P)  may  be  obtained 
directly  and/or  by  solving  other  related  problem(s);  it  depends  on  the 
definition  and  the  nature  of  the  components  comprising  these  data  elements 
for  a  specific  application  area. 


2.  DEVELOPMENT  OF  THE  SOLUTION  ALGORITHM 


The  solution  algorithm  that  has  been  developed  ti  solve  problem  (P) 
is  a  branch-and-bound  procedure  which  makes  use  of  Lagrangian  relaxation 
as  a  basic  step. 

This  chapter  considers  two  different  Lagrangian  relaxations  of 
problem  (P),  their  general  characteristics,  and  some  useful  results  leading 
to  the  specific  case  of  Lagrangian  relaxation  utilized  in  the  solution 
algor  itlun. 

2 . 1  Lagrangian  Relaxation 

Taking  a  set  of  "complicating"  constraints  of  a  ganeral  mixed- 
integer  program  into  the  objective  function  in  a  Lagrangian  fashion  (with 
fixed  multipliers)  results  in  a  "Lagrangian  relaxation"  of  the  original 
problem  [Geoff rion  (1974)].  The  relaxed  problem  is  easy  to  solve  compared 
to  the  original  problem,  and  provides  a  lower  bound  (for  minimization 
problems)  on  the  optimal  value  of  the  original  problem. 

Although  the  use  of  Lagrangian  relaxation  in  discrete  optimization 
has  been  reported  prior  to  1970  [e.g.,  Lorie  and  Savage  (1955),  Everett 
(1963),  and  Gilmore  and  Gomory  (1963)],  the  "birth"  of  the  Lagrangian 
approach  as  it  exists  today  [Fisher  (1978)]  occurred  in  1970  with  the 
successful  application  of  Lagrangian  relaxations  to  the  traveling  salesman 
problem  [Held  and  Karp  (1970,  1971)].  This  was  followed  by  application 
of  Lagrangian  relaxation  to  scheduling  problems  [Fisher  and  Schrage  (1972),  and 
Fisher  (1973,  1976)],  the  general  integer  programming  problem  [Shapiro 
(1971),  and  Fisher  and  Shapiro  (1974)]  and  the  generalized  assignment  problem 
[Ross  and  Soland  (1975)].  Table  2  lists  the  applications  of  Lagrangian 
relaxation  as  given  by  Fisher  (1978).  A  review  of  Lagrangian  relaxation 
is  also  provided  by  Shapiro  (1977)  and  C'nr istof ides  (1980). 
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TABLE  2 

★ 

APPLICATIONS  OF  LACRANGIAN  RELAXATION 


Problem 

Researchers 

Lagrangjan  Problem 

TRAVELING  SALESMAN 

Symmetric 

Asymmetric 

Symmetric 

Asymmetric 

SCHEDULING 

n|m  Weighted 
Tardiness 

1  Machine  Weight 
Tardiness 

Power  Generation 
Systems 

GENERAL  IP 
Unbounded  Variables 
Unbounded  Variables 
0-1  Variables 

LOCATION 

Uncapacitated 

Capacitated 

Databases  in 

Computer  Networks 

GENERALIZED  ASSIGNMENT 


Held  &  Karp  (1970,  1971) 
Bazarra  &  Goode  (1977) 

Balas  &  Christofides  (1976) 
Balas  &  Christofides  (1976) 


Fisher  (1973) 

Fisher  (1976) 

Muckstadt  &  Koenig  (1977) 

Fisher  &  Shapiro  (1974) 
Burdet  &  Johnson  (1976) 
Etcheberry,  et.  al.  (1978) 

Cornuejols,  Fisher,  & 
Nemhauser  (1977) 

Geoffrion  &  McBride  (1977) 
Fisher  &  Kochbaum  (1978) 

Ross  &  Soland  (1975) 
Chalmet  &  Gelders  (1976) 


SET  COVERING— PARTITIONING 

Covering  Etcheberry  (1977) 

Partitioning  Nemhauser  6  Weber  (1978) 


Spanning  Tree 
Spanning  Tree 
Perfect  2-Matching 
Assignment 


Pseudo-Polynomial 
Dynamic  Programming 
Pseudo-Polynomial  DP 
Pseudo-Polynomial  DP 


Group  Problem 
Group  Problem 
0  -  1  GUB 

0-1  VUB 
0  -  1  VUB 

0  -  1  VUB 


Knapsack 

Knapsack,  0-1  GUB 

0  -  1  GUB 
Matching 


*Source:  Fisher  (1978) 


-  17 


2.1.1  Relaxing  Problem  (P) 

Bv  dividing  constraints  (5)  by  s  ami  letting  r.  =  d...  ,  , 

k  tjk  1 1  '  s  k 

problem  (P)  can  be  restated  as  follows. 


Minimize  Z  Z  a,,  x..  +  £  b,  y, 

i-1  j=i  U  k==a  k  k 


sub  ject  to  Z  x . .  =  1 
i=l 


J=- . n 


X  **  riik  Xii  -  yk  k=1 . . 

i=l  j=l  1JK  1J  k 


x..,  y,  =  0  or  1  for  all  i,j ,k 
ij  k 


A  Lagrangian  relaxation  (LR^)  problem  (P)  relative  to  constraints 


(2)  is  obtained  as 


Minimize  Z  Z  a.,  x..  + 

i=l  j=l  1J  1-,  k 


(LK*  subject  to  _Z^  _^r..k  x..  <  yR 


P  n  /  m  \ 

^  -  v 


x  .  . ,  y.  =  0  or  1 
ij  'k 


k=l, . . . ,p 


for  all  i  ,  j  ,  k 


wliere  the  u  are  Lagrange  multipliers;  it  follows  that  the 

optimal  value  of  problem  (LR^)  is  a  lower  bound  on  the  optimal  value  of 

problem  (P),  i.e.,  Z(LR  )  <  Z(P)  .  We  will  continue  to  use  this  notation 

u  — 


in  which  /.(•)  is  the  optimal  value  of  problem  (•)• 


Another  I.agrangian  relaxation  (LRv)  of  problem  (P),  relative  to 


constraints  (5'),  is  obtained  as 
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Minin'  ize 


E  E  a  x .  .  +  E  b  y  -  E 
i]  ‘J  11  k  k  k  k 


\  K  ■  \  5  xu  rijk) 


subject  to  (2)  and  (6),  or  equivalently. 


(LR  ) 


Minimize 


Z  E  x  .  ( 

i  j  1J  \ 


a .  .  +  E  v,  r  .  . 
U  k  R  H 


ik)  -  J  y:<  (\  -  \) 


subj  ect  to  E  x . .  =  1 


i  1J 


j=3 . p 


x..,  y,  =  0  or  1  for  all  i,j,k 
ij  k 


(13) 

(2) 

(6) 


where  the  are  non-negative  Lagrange  multipliers;  it  follows 

that  Z(LR  )  <  Z(P)  , 
v  — 


2.1.2  General  Characteristics 

A  Lagrengian  relaxation  provides  a  lower  bound  on  the  optimal  value 
of  the  orginal  problem,  i.e.,  in  our  case  Z(LR^)  <  Z(P)  and  Z(LR^)  <_  Z(P) 
The  usefulness  of  a  Lagrangian  relaxation  depends  on  the  closeness  of  this 
lower  bound  to  the  optimal  value  of  the  orginal  problem.  However,  the 
relaxation  must  be  "easy"  to  solve  relative  to  the  original  problem.  We 
observe  that  the  optimal  value  of  y^  in  problem  (LR^  is  1  if 

(v,  -  b, )  >  0  and  0  if  (v,  -  b,  )  <  0  ,  and  then  prol lem  (LR  )  reduces 
k  k  —  k  k  —  v 

to  n  0-1  "multiple  choice"  problems  which  are  very  easy  to  solve. 

On  the  other  hand,  problem  (LR  )  reduces  to  k  0-1  knapsack  problems. 

u 

However,  these  problems  are  not  independent  because  of  the  interaction 
of  constraints  (5')  and  the  indirect  relationship  described  earlier  in 
Section  1.3  between  the  assignments  and  the  facilities.  In  view  of  this 
complexity,  relaxation  (LR^)  will  not  be  considered  further. 

The  choice  of  Lagrange  multipliers  in  relaxation  (LRv>  should 

be  such  that  ZO.R^)  is  as  large  as  possible  and  hence  as  close  as 

possible  to  Z(P)  in  view  of  the  relationship  Z(Lky)  _<  Z(P)  .  In  other 

words,  an  equivalent  problem  is  to  find  a  vector  v  (representing  v^,  v 

...,vk)  to 
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(D) 


Maximize 

v  >  0 


[Z(IV] 


(14) 


Obviously,  Z(LR  )  <  Z(D)  <  Z(P) 


The  general  properties  of  Lagrangian  relaxation  have  been  well 
described  in  the  literature  [e.g.,  Geoff rion  (1974),  Ceoffrion  and  McBride 
(197b),  and  Fisher  (1980)],  Some  of  these  properties  relating  the 
Lagrangian  relaxation  and  the  usual  LP  relaxation  are  stated  below. 

The  LP  relaxation  (P)  of  problem  (P)  is  obtained  by  relaxing  the 


integrality  constraints 

(6),  i.e. 

,  the  formulation  (P)  is 

r 

Minimize 

m 

E 

i=l 

n 

Z  a  .  x 

j-i  1J 

P 

.  .  +  X  by 

U  k=1  k  k 

(4) 

m 

subject  to 

Z  x  = 
i=i  iJ 

1  j=l,...,n 

(2) 

(!’)  4 

m 

n 

E 

i=l 

Z  r. 
j-i  ijk 

Xij  -  yk  k=1 . . 

(5’) 

k=l, .  . ,p 

(15) 

xij’ 

’k-1 

for  all  i , j , k 

(16) 

Note  that  the 

constraints 

x„  <  1  are  implicit  in  constraints  (2) 

Also  consider 

the 

following  partial  convex  hull  relaxation 

<P*> 

of  probLem  (P) . 

f 

m 

n 

P 

Mini mize 

E 

i=l 

Z  a  .  .  x 
i=i  1J 

.  .  +  E  by 

u  k-1  k  k 

(4) 

(P*>  - 

subject  to 

m 

E 

i=l 

n 

Z  r  .  . ,  x 
j=i  1)k 

Ij  i  yk  k'1 . . 

(5’) 

x .  . 
ij 

,  y^  £  convex  hull  {(2),  (6)} 

(17) 
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Then  tin*  relationships  between  the  optimal  values  of  various 
problems  [Ceuffrion  and  McBride  (1978)]  are  as  follows. 

Z(P)  <  Z(LR^)  <  max  Z(LR  )  =  Z(D)  *=  Z(P*)  <  Z(P)  (18) 

v  >  0 

a.  %  a,  % 

where  v  are  the  values  v^,  V2»...,Vp  of  a  dual  optimal  solution 
of  (P)  corresponding  to  constraints  (5'). 

Thus,  the  optimal  dual  solution  associated  with  the  usual  LP  relax¬ 
ation  furnishes  a  choice  of  Lagrange  multipliers  such  that  the  associated 
Lagrangian  relaxation  is  at  least  as  tight  as  the  usual  LP  relaxation, 
and  generally  a  good  deal  tighter  and  even  as  tight  as  the  partial  convex 
hull  relaxation. 

Since  Z(D)  =  Z(P*)  ,  the  quality  of  the  bound  obtained  from  the 

Lagrangian  relaxation  depends  on  where  Z(P*)  lies  in  the  range  between 
Z(P)  and  Z(P)  .  It  turns  out  that  problem  (LR  )  possesses  the 
"integrality  property,"  i.e.,  the  optimal  value  of  problem  (LK^)  is  not 

altered  by  dropping  the  integrality  conditions  on  its  variables  and 
therefore  (Ceoffrion  (197 A ) ] 

Z (D)  =  Z(P*)  =  Z(P)  (19) 

Thus,  the  Lagrangian  relaxation  (LRy)  is  no  better  than  the  LP 

relaxation  (P) .  On  the  other  hand,  Lagrangian  relaxation  (LR  )  does 

not  possess  the  integrality  property  and,  hence,  could  provide  an  equal 
or  better  bound  than  the  LP  relaxation  (P)  ;  but  the  computational 
difficulties  do  not  favor  pursuing  formulation  (LR^)  . 

It  is  possible  to  consider  alternative  formulations  of  problem 
(P)  with,  the  objective  of  obtaining  tighter  bounds.  This  aspect  is 
discussed  in  Chapter  6. 
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2.2  Sumo  lie_su_lj.  s 

We  now  turn  to  the  basic  question  of  choosing  Lagrange  multipliers 

v  so  that  (I.K  )  is  optimal  to  the  extent  possible,  which  is  equivalent 

to  solving  problem  (D) .  We  also  need  to  consider  this  question  when  some 

of  the  x..  and  y.  variables  have  been  assigned  values  of  1  or  0,  i.e., 
ij  k 

at  a  node  other  t can  the  starting  or  "root"  node  in  the  branch-and-bound 
tree.  For  this  purpose,  some  terminology  is  defined  and  formulations 
corresponding  to  problems  (P),  (LK  )  and  (D)  are  first  developed.  Then 

some  important  results  pertaining  to  the  choice  of  Lagrange  multipliers 
will  be  proved.  (lavish  (1978)  provides  a  method  of  obtaining  the  'best' 
multipliers,  based  on  solving  an  equivalent  linear  programming  problem. 
Such  a  formulation  is  difficult  in  our  case,  and,  besides,  we  propose 
to  avoid  solving  LP  problems  in  our  branch-and-bound  procedure. 

Define  the  sets 


s  =  { (i,j) i has  an  assigned  value  of  1  or  0}  ,  and 
T  =  {ic|y^  has  an  assigned  value  of  1  or  0)  . 


These  sets  represent  the  partial  solution  of  problem  (?)  and  the  variables 
contained  in  these  sets  are  termed  f ixed  variables.  [Geoff rion  (1967)]. 
Let  S  anti  T  represent  the  corresponding  complementary  sets,  i.e.. 


comprised  of  the  x. 


ij 


and 


y^  variables,  which  have  not  been  assigned 


specific  values  ;  ml,  tiierefore,  are  called  free  variables.  A  completion 

of  a  partial  solution  is  defined  as  a  solution  that  is  determined  by 

S  and  T  together  with  a  binary  specification  (0  or  1)  of  the  values 

oi  the  free  x.  and  y,  variables  from  sets  S  and  T  . 
ii  'L 


Let  SUS  =  Sj  and  TIJT  =  T  . 

Consider  a  partial  solution  to  problem  (P)  in  which  specific  values 
(of  i  or  0)  are  assigned  to  some  of  the 


x..  and  y  such  that 
ij  3  k 
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m 

Z  x  <  1  Vj  , 

i=l  J 
(i,j)eS 

and  E  X  r  x..  <  y,  VkeT 

ilk  xi  —  k 

1  j 

(i,j)cS 

E  Z  r .  x . .  <  1  Vkef 

i  j  ljk  ij  " 

(i, j)eS 


anu  sue 


:h  that  x.  =  1  ar.d  e  =  1  imply  tnat  keT  and  y,  =  1. 


Recall  that,  by  definition,  eik  -  1  if  design  i  uses  facility  k  , 

and  e  ,  =  0  otherwise, 
lk 

The  problem  of  finding  an  optimal  completion  of  the  partial 
solution  of  problem  (P)  can  be  stated  as  follows. 


Minimize 


l  l  a.,  *.,  +  r  b  y  +n  a,  X  4  +  Z  b.y 


i  J 

(i,j)eS 


ij  ij 


l  J  ij  ij  k  k  k  (20) 


keT  (i, j )eS 


subject  to 


Z  x ,  .  =  1  -  Z  x 
i  o  j  u 

(i,j)es  U,j)ES 


l  J  rljk  *ij  -  yk  "  *  j  cijk  *ij  Vk 


(i,j)eS 


(i,j)es 


x, . ,  y,  =  0  or  1  V(i,j)eS  ,  keT 

ij  k 


We  call  this  problem  (P^)  where  H  indicates  the  node  in  the  branch-and- 
bound  tree. 


A  Lagrangian  relaxation  of  problem  (P^)  with  respect  to  constraint: 
(22)  is  obtained  hy  introducing  non-negative  Lagrange  multipliers  , 


k*l,2,...,p  ;  the  relaxation  is  then 
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Minimi ze 


Z  Y  a  .  .  x  .  .  +  I  fa.  y.  i  H  a .  .  x .  .  +  Z  b,  y, 
ij.1J  1J  k  k  k  i  j  ^  k  k  k 


(i,j)eS 


keT  (i,j)£i 


Zv  y  -EL  r .  x . .  -  Z  Z  r  x .  . 

■  k  k  .  .  lj  k  ij  .  .  1  ]  k  ij 

k  1  j  i-  j 

(i,j)£S  (i, j)eS 


sub  j  ec  t  to  Z  x . .  =  1  -  Z  x  , 

i  i  J 

(i.jHS  (i,  j)eS 


x  .  .  ,  y,  =  0  or  1 
lj  k 


V(i,j)eS  ,  keT 


Rearranging  (24),  and  using  the  relationship  T  =  TUT,  we  have  problem 


V  : 


I  Minimize 


UR  „) 


Z  Z  x  /a  .  +  Z  v,  r  .  \  +  Z  Z  x.  .  ( a. .  +  Z  v,  r .  | 

1  J  (  1J  keT  k  Jk  I  i  j  Ul  U  kET  k  ^ 
(i, j)£S  '  1  '  (i, j)cs  '  1 

Vk  (Vk  ’  bk)  ^  (Vk  ■  tk)  (2M 


subject  to  l.  x..  =  1  -  Z  x.. 


(i, j)cS 


(i, j)es 


x .  .  ,  y  =  0  or  1 
ij  'k 


V(i,j)eS  ,  keT 


Then  we  have  Z(LR^  Z ( P ^  >  .  An  important  problem  is  the  choice  of 


Lagrange  multipliers  v^,v^,...,  ,  represented  by  vector  v  ,  that 

maximize  Z(LK^  y)  ,  i.e.,  the  problem  (D^)  : 


Maximize  £z(LK,  )j 


v  '•  0 
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We  now  state  and  prove  some  theorems  related  to  the  choice  of  Lagrange 
multipliers  v^»v^,...,v  . 


Theorem  1:  There  exists  an  optimal  solut  .on  to  problem  (D)  in 
which  -or  k  • 


Proof : 


Suppose  ,  in  an  optimal  solution  to  problem 

(D),  i.e.,  Z (D)  =  Z(LR  .)  where  /  *  <  b.  . 

V*  1  X 


Recall  that 

Z 


(LR  .  )  =  Min  Z  Z  x .  .  ( a .  .  +  Z  v*  r . \  -  Z  y,  (v*  -  b.  ^ 

v*  i  j  v  v 13  k  k  ljk;  k  k  [ k  v 


s.  t .  Z  x .  .  =  1 

i  iJ 


Vi 


Xij’  yk  "  0  or  1  Vi.j.k 


(2) 

(6) 


For  v*  <  b^  ,  the  optimal  value  of  is  0  ,  and  the  term  -  y^  (v*  -  b^ 

in  the  objective  function  is  0. 


Consider  what  happens  if  we  increase  v*  to  .  Call  the 

resulting  vector  v^  .  Consider  problem  (LRv)  •  The  optimal  value  of 

y^  in  problem  (LR^)  is  0  or  1  ,  and  the  term  -  y  1  _  ®  * 

However,  the  optimal  value  of  y^  is  the  same  in  problems  (LR^*)  and 

(LRv)  for  all  k  >  1  .  Therefore,  the  quantity  Z  y^  (v^  -  b^)  is  the 
—  k 

same  at  the  optimal  solution  for  both  v  =  v*  and  v  =  v  . 


Since  v^  >  v*  ,  we  note  that  in  the  objective  function. 


P  P 

a  +  Z  v  r  >  a.,  +  Z  v*  r 
ij  k=1  -k  ijk  -  lj  k  ijk 


Vi, j  , 


k=l 


and  therefore  Z(LRv>  Z(LRV*) 
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It  tallows  that  there  is  an  optimal  solution  to  problem  (D) 
in  which  v^  >  . 

Since  the  choice  of  k=l  was  arbitrary,  the  same  result  holds  for 

any  value  of  k,  k=l,...,p  ;  hence,  there  exists  an  optimal  solution  to 

problem  (D)  in  which  v,  >  b  for  all  k  . 

k  —  k 

Theorem _ 2:  There  exists  an  optimal  solution  to  problem  (D,.)  in 

which  >  b^  if  (i)  keT  or  (ii)  keT  and  y^  =  0 


Proof : 


Suppose  v^  <  in  an  optimal  solution  to  problem 


(Dj )  ,  i.e.,  Z(D^)  =  Z(LR^  ^*)  where  v*  <  b^ 

Then  k=l  can  he  such  that  keT  or  keT  . 


Case  (i):  Let  keT  . 
Recall  that 


Z(LR0  *)  =  Min  Z  Z  x..  (a  .  +  Z  vf  r... 
*"V  i  j  _  1J  \  1J  k  k  ljk, 

(i, j)cS 


+  Z  Z  x .  .  (i 

ij  V 


a,  .  +  Z  v*  r .  . 


(i, j )eS 


“  E-  yk 
keT 


vk  "  bk 


k  ijk 


'  Z  yl 
kET 


s.t.  Z  x .  .  =  1  -  Z  x. . 

i  iJ  1  iJ 

(i,j)cS  (i,j)cS 

x .  . ,  y,  =  0  or  1 
ij  7k 


Vi  (21) 

V(i,j)ES  ,  kef  (23) 


For  v*  <  b^  ,  ar.d  keT  ,  the  optimal  value  of  y^ 
term  -  y^  (v*  -  b^)  in  the  objective  function  is  0  . 


is  0  and  the 


2  b 


Lot  v*  bo  increased  to  b  ;  call  the  resulting  vector  v 


Consider  problem  (LR.  )  .  The  optimal  value  of  y,  in  (LR.  )  is 

Jo ,  v  r  J i  Z,  v 

0  or  1  ,  then  the  term  -  y  (v^  -  b1>  is  0  .  For  k  >  1  ,  the 


optimal  value  of  y  being  the  same  in  (LR,  *)  and  (LR.  )  ,  we 
k  x.,v  £,v 


find  that  Z  y  (v  -  b  )  is  the  same  at  the  oi timal  solution  for  both 
keT  K  K  k 


v  =  vA  and  v  =  w_  .  But  v^  >  v^  ;  therefore 


3ij  +*  ^  rijk  1  3ij  ^"ijk  V(i>j)eS  and  (i,j)£S 


Hence,  Z(LR  )  >  Z(LR.  *)  ,  wherefrom  it  follows  that  there  is  an 
x.,_v  —  &,v 


optimal  solution  to  (D.)  in  which  v  >  b  .  S'.nce  the  choice  of 

JL  1 


k=l  was  arbitrary,  the  same  results  hold  for  any  value  of  k,  kef  . 


Hence,  there  exists  an  optimal  solution  to  problem  (0^ )  in  which 


v.  >  b  for  all  keT  . 
k  —  k 


Case  (ii):  Let  keT  and  y,  =  0 


Considering  problem  (LR,  *)  ,  for  k  =  1  ,  /*  <  b  ,  and 

V  1  1 


y^  =  0  ,  the  term  -  y^(v"  -  b^)  in  the  objective  function  is  0  . 

Increase  v^  to  b^  and  call  the  resulting  vector  v^  .  The 
term  -  y^(v^  ~  b^)  is  0  .  For  k  >  1  ,  the  optimal  values  of  y^ 


are  the  same  in  problems  (LR.  *)  and  (LR.  )  .  Therefore  Z 

*’V  kcT 


y^(v^  -  b^)  is  the  same  at  the  optimal  solution  for  both  v  =  v* 


1 


and  v  =  v  .  Since  v^  >  v*  , 


*u  \  ruk  -  “u  \L  t  rijk  V(M),:S  and  a'j)E§ 
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Therefore  Z(I.R,  )  >  Z(LR.  *)  .  It  follows  that  there  exists  an  optimal 
i ,  v  -  l ,  v 

solution  to  (D^)  i-n  which  v,  >  .  The  choice  of  k=l  being  arbitrary, 

the  same  results  hold  for  any  value  of  k  ,  keT  and  =  0  ;  which  proves 

case  (ii)  of  the  Theorem. 

It  may  be  added  that  there  is  another  possibility  which  complements 
case  (ii)  of  Theorem  2,  i.e.,  if  kcT  and  y  =  1  .  We  treat  this  possi¬ 
bility  as  a  conjecture  since  a  result  similar  to  the  one  above  could  not 
be  proved,  as  discussed  now. 


With  yx  =  1  and  v*  h  j  ,  we  observe  from  problem  (LR.  v*) 
that  for  a  solution  vector  X-'  (with  elements  x^.*)  an(i  i  (with 
elements  y^*,  .  . .  ,v*  j  y^*=  1  and  y.,* . y  *  =  0  or  1)  , 


Z(LR,  *)  =  E  E  x..*  a..+v*r...+Z  v*r...\ 

i  j  _1J  V  1J  1J  k>l  k 

(i>j)eS 

+  S  E  x .  .  *  ( a .  .  +  v  *  r  .  +  E  v ,*  r  .  . ,  \ 

i  j  ij  \  ij  1  ijl  k>i  k  ijk J 

(i,j)eS 

v  (VJ ' bk) ■  yi(v • bi)  -  y**(v"b0 


Since  v^*  <  b^  and  y^  =  i  the  term  -  y^(v^*  “  b^)  is  positive. 
If  we  raise  v^*  to  b^  ,  say  v^  ,  the  term  -  y^(v^  ~  b^)  is  0  . 

The  difference  between  Z(LR  *)  and  the  objective  function 


problem  (LR^ 

.v> 

with 

X=X  and 

Y=Y  is 

E  E  x..*  v* 

r  ... 

+  (b, 

-  v  *)  -  >: 

Ex  *  b 

i  i  1J  J 

ijl 

1 

1 

i 

j  ij 

(bl  -  - 

e  r, 

i  j 

x  .  .  * 
1J 

(bi  -  v> 

riji  . 

This  difference  can  be  either  negative  or  positive,  and  s o  we  cannot 
conclude  that  there  is  an  optimal  solution  to  problem  (D,-,)  in  which 


.  ^  «  *- W  ►*  49 
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.  Wo  boliovo  Lliis  conclusion  to  be  false. 

k  -k 

Theorem  1:  Let  (X  ,  Y  )  solve  problem  (LR  )  for  v,  =  b, 
- * —  v  k  k 

*  * 

for  all  k  .  If  (X  ,  Y  )  is  feasible  for  problem 

(P) ,  there  exists  an  optimal  solution  to  problem  (D) 

in  which  v,  =  b,  for  all  k  . 
k  k 

Proof :  In  view  of  Theorem  1,  there  exisLs  an  optimal  solution 

to  (D)  in  which,  v^  ^  for  all  k  ,  i.e.,  v  >  b  . 

Let  v  be  such  an  optimal  v  .  We  will  show  that 
Z(lRv)  Z(LR^)  ,  from  which  it  follows  that  v  =  b 

solves  problem  (D)  . 


Recall  that 


Z (LR  )  =  Min  Z  Z  x. . 


X,Y  i  j 


ij 


(au  +  rk  ^  -  \) 


s.  t,  I  X  .  .  =1 


IJ 


x . 

ij 


^k 


=  0  or  1 


Vj 

Vi, j ,  k 


(2) 

(6) 


Since  v  >  b  ,  y,  =  1  V k  is  an  optimal  choice. 
-  k 


Hence,  Z(LR  )  =  Min  Z 

V 

~  X  i 


E  E  x.,  /a..  +  Z  v  r  \  -  E  /  v  —  b,  J 

ij  1J  \ij  k-*  ^7  k  \  k  7 


s.  t . 

Z  X  .  =  1 
i  ^ 

VJ 

(2) 

x  =  0  or  1 

ij 

Vi,j 

(6a) 

Now  consider  (LR^)  ■ 

Since  v  =  b  , 

the  last  term  of  the  objective 

function  drops  out,  and  we  have 


Z ( LR,  )  =  Mi 
b  X 


in  X  X  x  (a  +  X  b,  r  j 
,Y  i  j  ^  \  iJ  k  k  ljk/ 


subject  to  (2)  and  (6) 


-  2b  - 

-  Min  £  £  x.. 

I  a  .  +  £  b  r 

X  i  j  ^ 

^  il  k  k  i jk 

sub j  ect 

to  (2)  and  (6a) 

=  £  £  x  .  .  *  fa 

+  £  b  r  i 

ij  M 

ij  k  k  U  k/ 

where  Xx  with,  elements  x .  .  *  is  the  minimizing  solution  vector 

ij 

which  satisfies  (2)  and  (6a)  . 


*  /c 

Now  (X  ,  Y  )  feasible  for  (P)  implies  that 


subject  to  (2)  and  (6a) 


subject  to  (2)  and  (6a) 


_  £  £  x .  .  *  la  .  .  +  I  r  b,  \ 

i  j  U  [  U  k  Hk  kj 

►  rk  “A  i  ij  'jk  k  rk  kj 

<  ££x..*  /  a .  .  +  £  r  b,  \  =  Z(LR,) 

'  ij  U  \1J  k  ljk  k/  % 


.  4.-J  - 
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by  (27),  or  Z(LRv)  ^  Z(LR^)  ;  it  follows  that  v  =  b  solves 
problem  (D) . 


2 . 3  Relaxation  (PR^) 

Theorems  1  and  3  are  useful  in  providing  a  choice  of  Lagrange 

multipliers  as  a  starting  point  in  solving  a  relaxation  of  problem  (P) 

at  the  root  node.  Theorem  2,  similar  to  Theorem  1.  provides  results  for 

a  partial  solution  of  problem  (P) ,  i.e.,  at  a  node  other  than  the  root 

node  where  some  of  the  x..  and  y,  have  been  fixed  at  1  or  0  . 

lj  k 

Theorem  1  is  important  in  pointing  out  that  a  certain  set  of  Lagrange 
multipliers  v  such  that  v.  >  b,  for  all  k  would  provide  an  optimal 

K  K. 

choice.  Theorem  3  narrows  this  choice  to  v  =  b  for  all  k  for  a 

k  c 

specific  situation,  i.e.,  when  the  resulting  solution  is  feasible  for 
problem  (?) . 

Letting  v^  =  for  all  k  ,  problem  (LRv)  becomes: 

f 


Minimize 

£  £  c  .  . 
i  j  « 

X  .  . 

1J 

(28) 

subject  to 

\  *ij  =  1 

V3 

(2) 

x.  .  =  0 
ij 

or  1 

Vi,j 

(6a) 

where  c  .  . 

ij 

=  a  .  .  +  £ 

b.  r  .  ..  . 

x  13k 

(29) 

Note  that  problem  (LR^)  is  very 
the  sum  of  the  minimum  (over  i) 


easy  to  solve;  its  optimal  value  is  just 
c^.  for  all  j  ,  i.e. , 


Z(LR  )  =  £  min  {c . , } 
D  i  i  1J 


(30) 


We  solve  this  problem  as  a  starting  point  at  the  root  node  in  our  branch- 
and-bound  procedure. 
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5ki  Max  ( ak«.  ,  ^ 


so  y  is  fixed  at  1  and  keT  if  6^  =  1 


To  return  to  the  relaxation,  problem  (LR^  v) ,  we  must  make  a  choice 

of  the  vector  v  of  Lagrange  multipliers.  Of  course,  we  would  like  to 
use  an  optimal  choice,  i.e.,  a  vector  v  that  solves  problem  (D^)  • 

Recall,  however,  that  Theorem  2  did  not  provide  us  any  useful  information 
about  the  optimal  value  of  v^  if  keT  and  y^  “  1  .  To  simplify  our 

approacn  and  have  recourse  to  the  results  of  Theorems  1  and  3,  we  choose 

v  =  0  if  keT  and  y  =  1  .  Note  that  there  are  no  keT  such  that 
k  k 

y  =  0  because  of  practical  considerations  and  because  our  branching  rule 
only  results  in  fixing  y^  values  at  1  .  Problem  (LR^  v)  now  takes  the 


Minimise  I  I  x  .  .  (  a ,,+  Z  v,  r . . .  I  y,  /v,  -  b,  V  E 
i  j  1J  \  «  kef  k  ljkj  keT  k  (  k  V  keT 


V35) 


(..R^-)  i  subject  to  Z  x.j  =  1  , 


*ij,  yk  =  0  or  1  for  all  (i,j)e?  ,  keT  . 


Note  that  in  this  problem  (LR.  -) ,  v  *=  0  if  keT  .  Also  note  how  closely 

JC  j  V  K 

it  resembles  problem  (LR^),  the  relaxation  at  the  root  node.  As  in  that 

case,  we  would  like  the  lower  bound  Z(LR.  -)  to  be  as  large  as  possible, 

x ,  v 


i.e.,  we  seek  v  to 


Maximize  ^Z(LR^  -)J 


v  >  0 


Because  of  the  close  similarity  of  problems  (LR^  -)  and  (LRv), 

it  is  possible  to  obtain  results  about  problem  (Dj)  that  are  analogous 

to  those  obtained  about  problem  (D).  We  state  these  results  as  Theorems 
A  and  5.  Their  proofs  are  omitted  because  they  follow  precisely  the 


prool s  ol  Theorems  1  anil  i,  respectively,  and  their  validity  follows 

from  the  tact  that  problem  (LR  -)  is  essentially  the  same  as  problem 

■  ,  v 

(I.KV)  but  involves  only  the  free  variables. 

Theorem  4:  There  exists  an  optimal  solution  to  problem  (Dj) 

in  which  v  -•  b,  tor  all  kcT  . 
k.  —  k 

'k  k 

Theorem  5:  Let  (X  ,  Y  )  solve  problem  (LR„  -)  for 

X  ,  V 

_  k  k 

vk  =  b^  for  all  k£T  .  If  (X  ,  Y  )  satisfies 

(5')  for  all  keT  ,  there  exists  an  optimal  solution 

to  problem  (Dt)  in  which  v,  =  b  for  all  kef  . 
x,  k  k 

Just  as  Theorems  1  and  3  motivated  us  to  use  the  relaxation 
problem  (I.R^)  to  obtain  our  lower  bound  at  node  1  ,  Theorems  4  and  5 

motivate  us  to  set  for  ail  kef  in  relaxation  problem  (LR^  -) 

to  obtain  our  lower  bound  at  node  £  .  With  this  specification,  problem 
(LR  -)  becomes 

X,  V 


Minimize  !  Y.  a .  .  „  x  +  FC„ 

.  .  ijl  ij  £ 

i  J 


(PR, )  -I  subject  to  Y.  x.,  =  1 


x . .  =  0  or  1 
ij 


V (i, j)eS  , 


(23a) 


where 


c  .  .  „  =  a  .  .  +  Y  b.  r  .  . , 

lj£  1J  kef  k 


b,  ( 

k-i  k  v 


1  -  6 


k£ /  i  jk 


and  the  fixed  cost  FC^  is  given  by 


FC.  =  Z  b .  =  Z  6,  „  b.  . 
*  keT  k  k=l  kl  k 


(39) 
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This  specific  relaxation,  problem  (PR^) ,  is  01  the  same  form  as 
problem  (LR,  )  and  is  equally  easy  to  solve  in  one  pass.  Its  opLimal 
value  R(PR^)  serves  as  the  lower  bound  at  node  P  .  Note  that  for  £=1 
problem  (PR^)  is  the  same  as  problem  (LR^)  . 

It  is  clear  that  setting  each  Lagrange  multiplier  v^  to  b^ 

for  keT  and  to  0  for  keT  is  not  generally  optimal  in  terms  of 
achieving  the  tightest  lower  bound  (except  as  per  Theorem  3) .  But  it 
provides  a  good  starting  point  in  seeking  an  optimal  vector  v  and  it 
provides  an  easily  calculated  lower  bound  at  each  node  of  our  branch- 
and-bound  procedure.  The  question  of  how  to  improve  upon  this  choice 
of  multiplier  values  will  be  discussed  in  Chapter  6. 


3 .  MET’.'.ODOLOCY  FRAMEWORK 


The  branch-and-bound  procedure/methodol ogy  developed  to  solve 
problem  (P)  uses  Lagrangian  relaxation  (PR^)  as  a  basic  step.  The 

branching  rule  dictates  which  x..  variable  to  branch  on  at  each  node. 

In  addition,  there  are  certain  rules  (e.g.,  the  capacity  rule  and  the 
bounding  rule)  which  contribute,  significantly,  in  improving  the  overall 
efficiency  of  the  procedure. 


Some  basic  terms  such  as  fixed  and  free  variables,  partial  solution 
and  its  completion  were  introduced  in  the  previous  chapter.  This  chapter 
first  provides  a  preliminary  discussion  of  the  branch-and-bound  methodology, 
[Ceoffrion  (1967),  and  Geoff rion  and  Ifarsten  (1972)].  Representation  and 
storage  of  the  x_  variables  for  branching  and  backtracking  is  described 

in  order  to  pro\ ide  continuity  and  consistency  with  the  computer  program 
covered  in  Chapter  4.  This  is  followed  by  a  description  of  the  major 
components  of  the  branch-and-bound  methodology. 


Branching  and  backtracking  is 

branching  commences  by  fixing  the  x 

rule)  to  l  and  moving  to  the  left 
1  ix  the  corresponding  x..  variable 


done  on  the  x..  variables.  The 
ij 

variable  (selected  by  the  branching 

branch  node.  When  backtracking,  we 
at  0  and  move  to  the  right  branch 


node  (if  the  rignt  branch  node  has  not  already  been  explored).  An  x^j 

variable  can  also  bo  fixed  at  0  or  1  by  rules  other  than  the  brandling 
rule.  The  capacity  rule  and  the  bounding  rule  are  two  such  rules  employed 
in  our  methodology. 


Figure  4a  shows  a  branch-and-bound  tree.  The  x  variables  fixed 

at  0  or  1  at  any  node  due  to  rules  other  than  the  brandling  rule  arc 
shown  in  parenthesis  at  the  appropriate  node. 
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:<ode  ©  is  the  root  node  and  also  the  parent  node  for  nodes 
©  and  © 

Node  ©  is  the  parent  node  for  nodes  ©  and  ©,  etc. 


Figure  4a. 

A  branch-and-bound  tree  illustration 


Node  (ft) 

1 

2 

3 

4 

5 

6 
7 


Partial  Solution  (S^) 


<P 

{103 v  -  301,  -  401} 

{103,  -  301,  -  401,  402} 

(103,  -  301,  -  401,  -  402, 

{103,  -  301,  -  401,  -  402, 

{103,  -  301,  -  401,  -  402, 

{-  103} 


201,  101} 

201,  101,  204} 
201,  101 ,-204} 


Partial  solutions  for 


Figure  4b. 

the  above  illustration  (Figure  4a) 
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For  problem  (PR^)  ,  a  partial  solution  corresponding  to  set  S  at 

node  l  ,  i.e  ,  S,  contains  x..  variables  assigned  values  of  1  or 
i  xj 

0  .  For  simplicity  in  the  computer  program,  an  x^_,  variable  fixed  at 

1  is  represented  as  (100  i  +  j)  ,  whereas  an  x  variable  fixed  at 

0  as  -  (100  j  +  j)  ,  e.g.,  x  =  1  and  x32  =  0  a:e  represented  as 

302  and  -  302  respectively.  Since  branching  is  done  on  x  variables, 

it  is  necessary  to  make  a  distinction  between  x„  variables  fixed  at  1 

due  to  the  branching  rule  and  those  fixed  at  1  due  to  the  other  rules. 

We  make  this  distinction  by  underlining  the  positive  number  to  represent 

an  x  .  fixed  at  1  due  to  the  other  rules.  For  example,  204,  -  301, 
id 

_103  represent,  respectively,  x2^  =  1  due  to  the  branching  rule,  x^  =  0 
due  to  the  branching  rule  or  any  other  rule,  and  x^3  =  1  due  to  a  rule 
other  than  the  branching  rule. 


Figure  4b  shows  the  partial  solutions  ot  the  branch-and-bound 

tree  in  Figure  4a. 

lmpLic.it  enumeration  involves  generating  a  sequence  of  partial 
solutions  and  simultaneously  considering  all  completions  of  each.  For 
our  minimization  problem,  we  start  with  an  initial  solution  having  a 
very  large  value  (infinity)  as  an  initial  upper  bound.  As  the  computations 
proceed,  feasible  solutions  (those  satisfying  the  capacity  constraints) 
are  discovered  from  time  to  time,  and  the  best  one  yet  found  is  retained 
as  an  incumbent  solution  with  the  corresponding  value  as  the  best  upper 
bound .  It  may  happen  that  for  a  given  partial  solution  we  can 

determine  a  best  completion  of  ,  i.e.,  a  feasible  completion  that 

minimizes  the  objective  function  value  among  all  feasible  completions 
of  .  If  such  a  best  feasible  completion  is  better  than  the  best 

upper  bound,  then  it  replaces  the  latter.  Or  we  may  be  able  to  determine 
that  S  has  no  feasible  completion  better  than  the  incumbent.  In  either 

case,  we  can  f .1  thorn  .  (Various  situations  of  t ithoming  and  back- 
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tracking  in  our  branch-and-bound  procedure  are  described  in  the 
following  discussion.)  All  completions  of  a  fathomed  partial  solution 
have  been  implicitly  enumerated  in  the  sense  that  they  can  be 

excluded  from  further  consideration  (with  the  exception  of  the  relevant 
best  feasible  solution  of  if  it  has  been  retained  as  the  best 

upper  bound) . 

In  our  branch-and-bound  procedure,  at  any  given  node  where  we  can 
fathom  ,  we  backtrack  to  the  parent  node  and  move  to  the  right-hand 
branch  (if  that  branch  has  not  already  been  explored)  by  fixing  the 
appropriate  variable  at  0  .  However,  if  the  right-hand  branch 

has  already  been  explored,  we  continue  backtracking  to  a  parent  node  where 
we  can  move  to  a  right-hand  branch.  For  example,  in  Figure  4a,  when 
backtracking  from  node  3,  we  move  to  the  parent  rode  2,  and  to  the  right 
to  node  4  by  setting  x^  =  0  .  However,  when  backtracking  from  node  6, 

we  move  back  to  node  4,  then  back  to  node  2,  then  back  to  node  1,  and  to 
the  right  to  node  7  by  setting  x^  =  0  . 


On  the  other  hand,  if  the  partial  solution  S,  cannot  be  fathomed, 
we  branch  to  the  left  and  augment  S,,  by  fixing  a  free  variable  x..  at 

x,  1J 

1  (based  on  the  branching  rule),  and  then  we  try  to  fathom  the  resulting 

partial  solution.  In  addition  to  the  one  variable,  selected  by  the 

branching  rule,  some  other  free  x. .  variables  can  also  be  fixed  at 

ij 

0  or  1  according  to  the  application  of  rules  other  than  the  branching 
rule.  Note  that  this  can  also  happen  when  backtracking,  i.e.,  when 


has  been  fathomed  and  we  backtrack  and  move  to  the  right  by  setting 

the  appropriate  x..  variable  to  0  . 
ij 


Let  us  consider  examples  of  both  situations,  i.e.,  when 
has  not  been  fathomed  and  when  has  been  fathomed.  In  Figure  4a 

we  cannot  fathom  (i.e.,  S  at  node  1),  so  we  move  to  node  2  by 
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augmenting  S.  by  t  ixing  x^  =  1  based  on  the  branching  rule,  and 

by  fixing  x^  =  0  and  x^  =  0  based  on  the  application  of  the  other 

rules.  Similarly,  we  move  from  node  2  to  node  3  by  augmenting  by 

fixing  x  =  L  .  As  an  example  of  backtracking,  when  we  fathom  , 

we  move  back  to  the  parent  node  2,  and  to  the  right  to  node  4,  getting 
a  new  partial  solution  S,  by  replacing  x,  0  =  1  with  x. ,  =  0  ,  and 

further  augmenting  it  by  fixing  x.,.  =  0  and  x..  -  1  based  on  the 

application  of  the  other  rules. 

Computationally,  the  storage  and  update  of  partial  solution 
is  easily  accomplished  by  considering  Figure  4b.  If,  at  a  given  node, 
the  parLial  solution  has  not  been  fathomed,  e.g.,  at  node  4, 

determine  the  next  branching  variable  by  using  the  branching  rule, 
i.e.,  x.-,^  ,  and  augment  by  adding  204  as  the  last  entry.  Also, 

augment  S/(  with  any  other  free  variables,  if  appropriate, 

depending  on  the  application  of  the  other  rules.  Now,  consider  the 
case  where  the  partial  solution  has  been  fathomed,  e.g.,  at  node  6, 

and  we  backtrack;  starting  with  the  last  entry  in  Sf  .  we  consider  one 

entry  at  a  time,  going  backwards,  until  we  find  a  posi'ive  number  which 
is  not  underlined.  in  our  example,  it  is  103.  In  other  words,  we  must 
branch  to  tie  right  by  fixing  X  =  0  ,  i.e.,  we  replace  103  with  -  103 

and  we  are  at  node  7,  Should  we  find  that  we  have  no  positive  number, 
the  procedure  terminates  since  we  are  back  at  the  root  node  and  the 
right  branch  has  already  been  explored.  This  happens  when  backtracking 
from  node  7. 

In  t  lie  branch-and-bound  procedure  we  generate  a  sequence  of 
partial  solutions  as  we  move  from  one  node  to  another.  This  sequence  is 
nett-redundant  in  the  sense  that  no  completion  of  a  partial  solution 
ever  duplicates  a  completion  of  a  previous  partial  solution  that  has 


been  fathomed. 
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Since  one  of  the  x  values,  for  each  j  ,  must  be  1  ,  a  total  of 

ij 

(2m-l)  nodes  are  theoretically  possible  for  complete  enumeration. 
However,  most  of  the  solutions  may  be  infeasible  because  of  the  capacity 
constraints.  The  branch-and-bound  procedure,  thro-.gh  a  judicious  choice 
of  branching  variables,  and  elimination  of  certain  infeasible  and  non- 
optimal  assignments  through  various  rules,  turns  out  to  be  a  practical 
and  computationally  efficent  algorithm.  The  various  components  of  this 
procedure  are  described  next.  Detailed  procedural  steps  and  the  solution 
of  a  test  problem  will  be  covered  in  Chapter  4. 

3 . 1  Bounds 

3.1.1  Lower  Bound 

At  a  given  node  Z  in  the  branch-and-bound  tree,  a  lower  bound 
(LOWS)  is  obtained  by  solving  relaxed  problem  (PR^)  . 


LOWB  =  Z (PR. ) 

£ 


(40) 


Recall  that  problem  (PR^)  is  very  easy  to  solve  by  considering  the 

minimum  c  .  . .  over  those  j's  for  which  x..  is  nit  fixed  at  1  , 

ij 

i.e.,  jeW  ,  where  W  is  the  complement  of  W  defined  by  expression  (32), 


Z(PR  )  =  t  c...  +  E 

*  jew  jew 


*"  cijt  +  FCt 
(i,  j)es 


(4i) 


where  given  by  expression  (38),  l.e., 


Cij£  =  aij  +  “  bk  (1  -  \l>  rijk  ’ 


(38) 


and  the  fixed  cost  (FC^)  is  given  by  expression  (39),  i.e.. 


FS  ' F  V  \  • 

k 


(39) 


where  6^  is  given  by  expression  (34). 


1 


I 


Note  that  it  none  of  the  x  variables  is  l'axed  at  1  ,  as  is 

ij 

generally  the  ease  at  the  root  node,  then  all  6,  ,  =  0  ,  and,  therefore, 

k-L 

P 

FC  -  0  ,  and  e  -  a..  +  Z  b  r  .  Z(PR  )  is,  then,  simplv  the 
1  ljl  ij  k=1  k  ljk  1  ' 

middle  part  ol  expression  (41).  Vie  use  the  term  "generally"  because 
it  is  possible  that  the  capacity  rule  could  forca  certain  x 

ij 

variables  to  1  (or  0)  at  the  root  node,  prior  to  solving  the  relaxed 
problem  (PR^)  . 

1.1.2  Upper  Bound 

At  any  given  node  1  ,  let  X  =  {x_^}  represent  the  solution  of 
problem  (PR^)  .  it  this  solution  is  feasible  for  problem  (P),  i.e., 
if  X  satisfies  the  capacity  constraints  (5)  or  (5') 

l  dijk  xij  -  \  yk  vk  ,  (42) 

where  y,  =  i  if  Z  1  d..,  x  >  0  , 

k  A  ■  tjk  ij 

X..CX  (43) 

ij 

=  0  otherwise  , 


then  the  value  of  problem  (P)  corresponding  to  this  solution 
gives  an  upper  bound  (1’PB)  : 


Z  X  a .  .  x .  .  +  Z  b  y,  , 
•  •  ij  ij  .  k  k 

i  J  k 

x.  .ex 

ij 


where  y  is  def ined  by  (4J)  . 


3.1.3  I’est  Upper  Board 


A  current  Lowest  upper  bound  is  retained  as  the  lest  upper  bound 
(BUB),  the  corresponding  solution  X  representing  the  incumbent  solution. 
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The  branch-and-bound  procedure  is  initiated  by  assuming  a  very  large 
value  as  the  best  upper  bound,  and  is  replaced  by  better  (lower)  values 
as  the  procedure  continues. 

A  positive  fractional  value  c  can  be  specified  if  a  sub-optimal 
solution  is  acceptable.  For  example,  for  E  =  0.001  ,  the  resulting 
solution  value  is  guaranteed  to  be  within  0.1  percent  of  the  optimal 
solution  value.  When  £  is  non-zero,  the  adjusted  best  upper  bound  (BUBS) 
is  defined  as: 

BUBS  =  BUB/ (1  +  £)  .  (45) 

Obviously  when  £  =  0  ,  BUBS  =  BUB  . 


3.2  Facility  Usage  Rule 

This  rule  is  used  to  identify  facilities  forced  into  usage  at  a 
given  node  l  and  hence  fix  corresponding  free  variables  y^  at  1  . 


For  a  partial  solution  ,  define 


V  =  dijk  lf  j£W  ’ 


=  min  d ,  . ,  if  jeVJ  . 
(i, j)e§ 


(46) 


The  facility  usage  rule  states 
is  not  already  fixed  at  1  ,  if 


that 
E  d 
j 


j  k£ 


for  any  facility  k  , 
>  0  ,  tnen  facility 


where 
k  is 


forced  into  usage  and,  therefore, 


should  be  fixed  at 


1  . 


This  rule  is  applied  at  every  node  prior  to  applying  the  capacity 
rule.  In  other  words,  this  rule  is  applicable  to  capacitated  as  well  as 
uncapacitated  problems. 

3 . 3  Capacity  Rule 

This  rule  is  designed  to  "exclude"  infeasible  assignments  prior 
to  solving  the  relaxed  problem  (PR^)  •  This  is  done  by  exploiting  the 


-  4  j 


rei.it ionship  between  the  capacities  required  (d .  )  and  the  capacities 

1  j  k 

available  (a,  )  for  a  given  partial  solution  oi  problem  (P). 

\ 

The  capacity  rule  states  that  for  a  facility  k  and  an  activity 

j  , "exclude" a  free  x..  variable  (i.e.,  fix  it  at  0)  for  which 
J  ij 


(dijk  -  >  <*k  -  ?  v>  - 


(i, j)s  S 


C47 ) 


where  d.,  .  is  defined  by  expression  (46).  The  right-hand  side  of  this 
J  kJC 

inequality  (47),  when  positive,  represents  the  available  capacity  at 
facility  k  .  The  left-hand  side  shows,  for  a  given  j  ,  the  difference 
between  a  d...  corresponding  to  a  free  x..  variable  and  d . ,  „  . 

ljk  c  &  ij  j  kx. 

If,  for  a  specific  d...  ,  this  difference  is  more  than  the  available 

ljk 

capacity,  the  corresponding  free  x..  variable,  if  fixed  at  1  ,  would 

result  in  an  infeasible  solution.  Thus,  by  looking  ahead,  we  can  exclude 
such  a  free  variable  by  assigning  it  a  value  of  0  . 

Note  that  if  the  right-hand  side  of  expression  (47)  is  negative, 
then  any  completion  of  such  a  partial  solution  will  be  infeasible  and 
we  backtrack  in  our  branch-and-bound  procedure. 

The  capacity  rule  is  applied  to  all  the  facilities  by  considering 
one  facility  at  a  time.  The  cycle  of  examining  all  tht  facilities 
continues  until  no  more  assignments  can  be  excluded.  During  the  course 
of  application  of  this  rule,  if  all  but  one  of  the  free  x^  variables 

have  been  excluded  (fixed  at  0)  for  a  given  j  ,  then  that  particular 

x..  variable  is  fixed  at  1  because  of  constraints  (2),  i.e.,  each 
ij 

activity  j  must  be  assigned  to  one  and  only  one  design  i  .  The 
partial  solution  is  updated  accordingly  to  reflect  the  x  variables 

fixed  at  0  or  1  due  to  the  application  of  the  capacity  rule. 


The  capacity  constraints  for  an  uncapacitated  problem  are  not 

Hence,  the  capacity  rule  is  useful  only  for  capacitated  problems. 


active. 
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3 . 4  Branch ing  Rule 

This  rule  provides  the  choice  of  the  x^.  variables  on  which  to 

branch.  If  the  partial  solution  at  a  given  node  SI  is  not  fathomed,  we 
branch  further  by  fixing  a  free  x  variable  at  1  and  moving  to  the 

left  branch  node. 

According  to  the  branching  rule  the  choice  of  the  branching 

variable  depends  on  the  c  .  values  and  is  such  that  the  corresponding 

ij  x. 

x..  ,  if  perturbed,  has  the  maximum  impact  on  the  optimal  value  of 


problem  (PR^)  • 


For  a  given  j  ,  define  c  ,  the  minimum  permissible  c  • 

i^jit-  ijx. 


the 

second 

smallest 

permissible 

Cii£  ’  i-e” 

ijj^ 

=  min 
i 

Cij* 

for 

jeW  and 

(i, j)es 

(48) 

i2jC 

=  min 

i 

Cij«- 

for 

j  eW  and 

(i,j)eS 

(49) 

For  each  jeW  ,  define  D.0=c.  ..  -c.  ..  .  (50) 

Our  branching  rule  states  that  a  free  x^  variable  corresponding  to 

c.  such  that  D.0  is  maximized  over  all  j  ,  is  selected  as  the 

x^Jx,  j  x, 

next  branching  variable  and  assigned  a  value  of  1. 

3. 5  Bounding  Rule 

This  rule  is  designed  to  "exclude”  certain  non-optimal  assignments. 
These  assignments  cannot  lead  to  an  optimal  solution  as  we  branch  from 
one  node  to  the  next  left  branch  node. 


The  bounding  rule  states  that  a  free  variable  should  be 

excluded  (by  assigning  it  the  value  0)  for  which 


(c  .  .  .  -  <:  .  .  ,)  >  (BUBS  -  LOWB)  for  jcW  and  (i,j)eS 


V 


(51) 


ij'- 

I. 

where  c  ,  ,  BL'BS,  and  LOWB  are  given  by  expressions  (48)  ,  (45), 

v*- 

and  (40),  respectively. 


Thus,  by  looking  ahead,  we  exclude  those  assignments  which  will 
provide  lower  bounds  higher  than  BUBS. 

The  bounding  rule  is  applied  to  each  jeW  just  prior  to  selecting 

the  x  variable  for  branching  to  the  left. 
lj 


As  in 
in  excluding 

a  given  JCW  , 


the  case  of  the  capacity  rule,  if  the  bounding  rule  results 
(fixing  at  0)  all  but  one  of  the  free  x..  variables  for 

then  that  particular  x  .  variable  is  fixed  at  1  .  Also 


the  partial  solution  is  updated  accordingly  to  reflect  the  x.. 


variables 


fixed  at  0  or  1  due  to  the  application  of  the  bounding  rule. 


3 . h  Bae kLrncking  Rules 


If  a  partial  solution  at  a  given  node  has  beer,  fathomed,  we  back¬ 
track.  The  backtracking  rules  are  typical  of  a  brancli-and-bound  procedure. 
In  addition,  the  application  of  the  capacity  rule  and  the  bounding  rule 
can  lead  to  backtracking.  The  criteria  for  backtracking  include  the 
to  1  lowing. 


(a)  When  applying  the  capacity  rule,  if  the  available 

capacity  given  by  the  right-hand  side  of  inequality 
(47)  is  negative,  i.e.,  (s^  -  E  cd^)  <  ^  > 

then  backtrack.  ^ 


(b)  If  LOWB  ’■  BL'BS,  then  backtrack.  Otherwise  compute 
LiiB  it  the  solution  is  feasible  in  problem  (P)  . 

Then  update  BUB  and  BUBS  if  UPB  <  BUB;  and  back¬ 
track  if  LOWB  =  BUBS. 

(c)  If  further  brandling  is  not  possible,  then  backtrack. 

This  can  happen  due  to  the  capacity  rule,  the  bounding 
rule,  or  the  branching  rule  if  the  updated  partial 
solution  is  such  that  no  further  branching  is  possible, 
i.e.,  x..  variables  arc  fixed  at  1  for  all  j  ,  or 

equ  iva  lent  Ly  ,  W  =  <t> 


When  any  of  the  backtracking  criteria  apply,  we  backtrack  to 
the  parent  node  and  move  to  the  right  branch  node  (if  the  right  branch 
has  not  already  been  explored)  by  fixing  the  appropriate  x  variable 

at  0  .  If  the  right  branch  has  already  been  explored,  we  continue  back¬ 
tracking  to  a  parent  node  where  we  can  move  to  a  right  branch  node. 

The  branch-and-bound  procedure  terminates  when  we  backtrack  to  the 
root  node  and  find  that  the  right  branch  node  has  already  been  explored. 


\ 


f 


r 


4.  COMPUTATIONAL  STEPS  AND  THE  COMPUTER  PROGRAM 


A  computer  program  called  ZIPCAF  (an  acronym  for  Zero-one  Integer 
Program  for  multiac Liv ity  multifacility  Capacity-Constrained  Assignment 
Problems)  implementing  the  branch-and-bound  methodology  has  been  developed. 

Detailed  procedural  steps  and  guidelines  to  use  the  computer 
program  are  described  in  a  separate  document  [Chhabra  and  Soland  (1980)] 
titled  "Program  Description  and  User's  Guide  for  2.IPCAP — a  Zero-one 
Integer  Program  to  solve  multiactivity  multifacility  Capacity-constrained 
Assignment  Problems."  Specifically,  the  document  includes: 

.  Problem  formulation  (P)  and  potential  areas  of 
appl icat ion 

,  Overall  flow  diagram  and  detailed  procedural  steps 
for  the  computer  program 

.  Program  listing  and  dictionary  of  the  symbolic  names. 

The  listing  includes  extensive  use  of  comment  cards 
to  explain  various  computational  steps. 

.  User  information  including 

schematic  diagram  of  the  deck  structure, 

detailed  instructions  for  the  job  control  (JCL) 
cards,  program  parameter  card,  program  options 
card,  and  the  various  other  input  data  cards. 

.  Three  test  problems  to  demonstrate  the  use  of  the 
program.  The  display  includes  coded  input  and 
annotated  outputs  reflecting  the  use  of  selected 
prog  Lam  options. 

As  mentioned  ear  ier,  ZIPCAP  is  primarily  designed  for  capacitated 
problems.  However,  uncapac itated  problems  can  be  solved  as  a  special 
case,  and  this  is  demonstrated  by  including  an  uncapacitated  test 
problem. 
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Because  of  the  extensive  coverage  of  the  program  description  and 
user  guidelines  in  the  above  document,  this  chapter  provides  only  an  over¬ 
view  of  the  computer  program,  including  an  overall  flow  diagram,  and  a 
summary  of  the  program  options,  in  order  to  provide  continuity  in  this 
document.  In  addition,  a  step-by-step  description  01  a  test  problen  is 
presented  to  demonstrate  the  use  of  the  various  components  of  the  branch- 
and-bound  methodology.  The  computer  printout  showing  step-by-step  details 
is  obtained  by  use  of  one  of  the  program  options.  The  use  of  this  option 
to  display  detailed  steps  in  this  document,  in  fact,  complements  the  use 
of  the  various  options  demonstrated  in  the  other  document. 

4 . 1  The  Program 

Figure  5  presents  a  simplified  flow  diagram  of  the  branch-and- 
bound  procedure.  The  major  computational  steps  for  the  computer  program 
are  numbered  in  circles.  These  steps  are  essentially  based  on  the 
methodology  components  described  in  the  previous  chapter.  A  step-by-step 
description  has  been  included  in  the  other  document  [Chhabra  and  Soland 
(1980) ] . 

The  computer  program  ZIPCAP  is  written  in  FORTRAN  IV,  and  has 
been  developed  and  tested  on  an  IBM  3031  at  the  George  Washington 
University.  The  program,  comprising  about  480  lines  is  currently 
dimensioned  for  a  maximum  problem  size  of  35  designs  (ro, ,  35  activities 
(n)  and  30  facilities  (p) .  The  program  size  to  execute  a  problem  has 
two  components:  one,  due  to  the  program  itself,  comprising  173  K  bytes, 
and  the  other  dependent  or.  the  dimensions  of  the  arrays  given  by  the 
following  functional  relationship. 

f(m,a,p)  =  4[(p+4)mn  +  (m+5)p+9n]  bytes 

The  computer  program  listed  in  the  other  document  has  since  been 
further  improved.  The  basic  improvement  has  been  the  addition  of  the 
facility  usage  rule.  This  rule,  as  described  in  Chapter  3,  is  applied 
both  to  capacitated  and  uncapacitated  problems  just  before  the  applica¬ 
tion  of  the  capacity  rule.  For  completeness  of  this  document,  a  revised 
program  listing  is  included  in  Appendix  A.  It  may  !e  mentioned  that 
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Figure  5. 

Simplified  flow  diagram  for  the  branch-and-bound  procedure 
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the  revised  program  solves  the  test  problems  included  in  the  other  document 
more  efficiently  —  in  less  time  and  in  fewer  notes  (with  an  average 
reduction  in  nodes  of  31  percent).  The  improvement  in  efficiency  seems 
to  result  from  the  "multiplicative"  effect  of  the  various  rules.  Another 
improvement  made  is  that  the  computer  printout  always  displays  the  node 
number  (IBNOD)  at  which  the  best  upper  bound  changes  (improves)  and  the 
corresponding  values  of  the  best  upper  bound  (BUB)  and  the  adjusted  best 
upper  bound  (BIBS). 

ZIPCAP  provides  numerous  options  to  the  program  user.  These  options, 
described  in  the  other  document,  are  summarized  in  Table  3. 

Option  ICAPR,  the  capacity  rule,  is  automatically  skipped  by  the 
program  when  solving  an  uncapacitated  problem.  Option  ISTEP,  the 
intermediate  steps'  listing,  even  when  skipped,  provides  information  on 
the  total  number  of  nodes  explored.  A  summary  listing  provides  necessary 
information  tc  construct  the  branch-and-bound  tree,  whereas  a  detailed 
listing  of  the  intermediate  steps  is  useful  when  changing  or  debugging 
the  program. 

Option  EPS,  the  optimal/suboptimal  solution,  provides  the 
flexibility  of  obtaining  a  suboptimal  value  guaranteed  to  be  within 
a  specified  fraction  of  the  optimal  value.  The  resulting  solution  may 
he  suboptimal  but  could  provide  a  considerable  saving  in  terms  of 
exploring  fewer  nodes  in  comparison  to  those  necessary  for  obtaining 
an  optimal  solution. 

Option  ET,  by  providing  important  information  at  a  specified 
elapsed  time,  is  useful  in  a  situation  where  the  total  time  allocated 
to  solve  a  problem  may  not  be  sufficient  and  the  program  terminates 
before  verifying  an  optimal  solution.  The  information  provided  by 
this  option  includes  an  updated  partial  solution  showing  the  x^_ 

variables  fixed  at  0  or  i,  at  the  current  node  being  explored  at 
the  specified  tire  ET.  By  looking  at  the  first  few  variables  displayed 
in  the  partial  solution  of  the  current  node,  it  is  possible  to  assess 
the  extent  ot  the  branch-and-bound  tree  explored  until  time  ET .  For 


■ . 


SUMMARY  OF  ZIPCAP  OPTIONS 


Skip  this  option 


example,  in  view  of  the  terminology  in  Figure  4b  (Chapter  3),  if,  at  an 

arbitrary  node,  the  first  term  of  the  partial  solution  is  positive, 

i.e.,  the  x  .  variable  has  value  1,  then  we  are  still  in  the  left 
ij 

half  of  the  total  branch-and -bound  tree.  If  the  first  term  is  negative, 
i.e.,  the  x  variable  has  value  0,  then  we  arc  in  the  right  half  of 

the  total  branch-and -bound  tree  and  have  explored  half  of  the  total 
(theoretical)  solutions  corresponding  to  the  left  half  of  the  tree. 

If  the  first  two  terms  are  negative,  i.e.,  the  fi'-st  two  x^ 

variables  have  value  0,  then  one  quarter  of  the  total  (theoretical) 
solutions  remain  to  be  explored,  since  we  are  in  the  next  right  half  of 
the  right  half  of  the  total  branch-and-bound  tree,  as  illustrated  in 
Figure  6. 


Each  i egment  represents  1/4 
of  the  total  (theoretical) 
solutions 

Each  segment  represents  1/2 
Left  half  of  Right  half  of  of  total  (theoretical) 

the  "total"  the  "total"  solutions 

branch-and-  branch-and- 

bound  tree  bound  tree 


Figure  6 

Illustration  for  estimating  the  extent  of  the  branch1 
and-bound  tree  explored 
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Recall  Iron'  Chapter  3,  that  a  total  of  (2m-l)n  nodes  are  theoretically 
possible.  Thus,  it  the  first  g  [g  <  (ra-l)n]  terms  at  an  arbitrary 

node  are  negative,  then  theoretically  about  [(2n~l)n/  g]  nodes  remain 
to  be  explored. 


4 , 2  An  Illustrative  Example 

We  consider  a  capacitated  test  problem  with  five  designs  (m) , 
four  activities  (n) ,  and  eight  facilities  (p)  to  demonstrate  the  use 
of  the  branch-and-bound  procedure  and  the  computer  program. 

The  computer  printout  for  this  problem  shoving  step-by-step 
details  for  a  couple  of  nodes  is  presented  in  Appendix  B. 

As  shown  in  the  beginning  of  the  printout,  the  options  selected 

are: 


.  1INPT  =  1,  i.e.,  list  the  input  data 

.  ICAPR  =  1,  i.e.,  use  the  capacity  rule 

.  ISTEP  =  2,  i.e.,  list  detailed  intermediate  steps 

.  IUNCAP  =  1,  since  this  is  a  capacitated  problem 

.  EPS  =0.0  implying  that  an  optimal  solution  is  desired 

.  ET  =  0.0  since  a  detailed  listing  of  intermediate  steps 

will  be  available. 


Following  the  listing  of  the  options,  input  data  listed  for  the  problem 
include  variable  costs  a..  ,  fixed  costs  b,  ,  available  capacities 

s^  ,  and  capacities  required  d  ^  values  are  generated 

by  the  computer  program. 


The  computer  program  follows  the  procedural  steps  marked  in  the 
flow  diagram  presented  in  Figure  5.  These  steps,  llong  with  the 
relevant  terminology  used  in  the  computer  printout,  are  described  below 
for  a  couple  of  nodes,  followed  by  a  complete  branch-ana -bound  tree  for 
this  problem.  As  mentioned  earlier,  a  dictionary  of  the  symbolic  names 
used  in  the  computer  program  is  included  in  the  other  document. 


*1 
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Node  1 

Step  1:  Initialize. 

Initialize  BUB  =  9999999.0,  and  since  EPS  =  0.0,  BUBS  =  BUB. 

Also  S  =  <p  and  W  =  ip  .  In  the  computer  printout,  vector 

FIX(J)  represents  the  set  VJ,  and  matrix  CX(I,J)  /epresents  both,  fixed 

and  fret'  x  .  variables.  In  the  CX(I,J)  matrix,  an  x.,  variable  fixed 
ij  ij 

at  1  or  0  is  represented  as  1  or  2,  respectively,  and  a  free  x  variable 

is  represented  by  the  value  0.  Initially,  all  the  x  variables  are  free 

as  shown  by  matrix  CX(I,J)  in  the  printout. 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l, 2 , . . . ,8  . 

In  the  printout,  itlND(J)  represents  d ^  ^  defined  by  expression  (46), 
and  MINSD  represents  E  d.y  .  As  shown  in  the  printoit,  MINSD  is  0  for 

j  J 

k=l,2,...,8  ,  and  so  the  facility  usage  rule  does  not  'orce  any  facilities 

into  usage;  and  as  shown  by  matrix  CX(I,J)  for  k=l,2,...8  ,  the  capacity 

rule  does  not  fax  any  x..  variables. 

-LJ 

Step  3:  Solve  the  relaxed  problem  (PR^)  . 

In  the  printout  FLB(K)  represents  6  ^  ,  given  by  expression  (34), 
for  computing  FC,  ,  and  C(I,J)  represents  c . . .  defined  by  expression 

l  J  X- 

(38).  Being  at  the  root  node,  i  =  1.  Further  the  solution  of  problem 
(PRp)  ,  i.e.,  X  =  {x..}  is  shown  in  the  printout  by  SOLX(J)  which 

X.  1J 

for  (I’R.)  is  X  =  {x,  ,  =  x ,  ~  =  x_..  =  x,  ,  =1}  . 

1  41  42  21  A  4 

Step  4:  Compute  the  lower  bound. 

The  expressions  (40)  and  (41),  i.e., 

L0WB  =  ZO’Rj,)  (4°) 

=  c  .  +  X  min  c  +  FC  . 

.  ljJ  .  -  .  ij?  *  (41) 

jcW  jr.W  i 


(i, j)eS 


-  55  - 


are  represented  in  the  printout  as 

LOWB  -  M1NSC  +  FC 

=  729839.3125  +  0  =  729839.3125 


Step  5:  Compare  LOWB  with  BUBS. 

Since  LOWB  <  BUBS,  go  to  Step  6 

Step  6:  Check  if  solution  X  is  feasible  in  problem  (P),  i.e., 
expression  (42)  is  satisfied. 


5  ?  dijt  -ij  t  vk 

X.  .ex 

ij 


Vk 


(42) 


In  the  printout,  NSUMD  represents  the  left-hand  side  of  this  inequality, 
and  for  each  k  ,  the  capacity  constraints  are  satisfied. 

Step  7:  Compute  the  upper  bound. 

UPB  is  given  by  expression  (44),  i.e., 


UPB  =  T.  Z  a  x  +  E  b  y  (44) 

i  j  J  1J  k  k  * 

x .  ex 

U 


In  the  printout,  the  corresponding  expression  is  represented  as 

UPB  =  NSUMA  +  FCUB 

=  678,502  +  101,000  =  779502.0  . 

Step  8:  Compare  UPB  with  BUB. 

Since  UPB  <  BUB,  go  to  Step  9. 

Step  9:  Set  BUB  =  779502.0  .  Since  EPS  =  0.0,  BUBS  =  BUB. 

Since  LOWB  <  BUBS,  go  to  Step  10. 

Step  10:  Left  branching  is  possible  since  W  =  <p  as  sho;/n  by 
vector  FIX(J);  go  to  Step  11. 

Step  11:  Apply  the  bounding  rule  and  the  branching  rule. 


ij 


According  to  our  bounding  rule,  a  free  x 


variable  is  excluded 
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(fixed  at  0)  for  which 

(c  ,  —  c .  . „ )  >  (BUBS  -  LOWB)  for  jeW  and  (i,j)eS  (51) 

lji.  i^Jl 

For  x  ,  (210,381.4375-145,201.5)  >  (779, 502 . 0-729, 839. 3125) . 

This  also  holds  for  x„_  and  x, ,  ,  i.e.,  the  bounding  rule  results 

33  14 

in  fixing  x^  >  x.^  >  and  x^  at  0.  This  is  shown  in  the  printout 

by  matrix  CX(I,J)  where  the  corresponding  variables  have  been  assigned 
the  value  2  because  of  the  bounding  rule. 

The  branching  rule  directs  us  to  select  a  free  x..  variable 
°  ij 

corresponding  to  c.  . „  for  which  D.s  =  c.  .  *  -  c.  .  p  is  maximized 

J  ~  J  ^  1*2  ~ 

over  all  j  .  In  the  printout,  c.  .  0  ,  c.  . «  ,  and  D . „  are  represented 

i^J^-  J*- 

by  NMINC(J),  MINC(J)  and  DIFBR(J),  respectively.  Since  D2^  is  the  maxi¬ 
mum,  is  selected  as  the  next  left  branching  variable.  This  is  shown 

in  the  printout  by  BR1  and  is  represented  as  (100  i  t  .' )  e.g.,  402  . 


Using  the  terminology  employed  in  Figures  4a  and  4b,  the  x_ 

variables  fixed  at  0  or  1  in  the  partial  solution  will  be  shown 

as  =  {-  103,  -  303,  -  104,  402}  .  In  the  computer  printout,  vector 

STX  displays  the  x..  variables  fixed  at  0  or  1.  The  representation 
ij 

of  the  variables  is,  however,  somewhat  different.  An  x..  variable 

13 

fixed  at  0,  due  to  any  rule,  is  shown  as  -  (100  i  +  j)  -  1,000,000, 
e.g.,  x,,  is  shown  as  -  1,000,103;  an  x..  variable  fixed  at  1  due 

to  the  brandling  rule  is  represented  as  (100  i  +  j),  e.g.,  x^  as  402; 

and  an  x..  variable  fixed  at  1  due  to  a  rule  other  than  the  branching 
ij 

rule  is  shown  as  (100  i  +  j)  +  1,000,000,  e.g.,  x^  is  represented  as 
1,000,203. 

In  the  printout,  vector  STX  represents  updated  partial  solution 


We  now  move  to  Node  2. 


Node  2 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
printout . 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l, 2, ... ,8  . 

As  shown  in  the  printout,  MINSD  (representing  Z  d  )  ,  being 

j  2 

positive  for  1=1,2, 3, A,  and  5  ,  these  facilities  are  forced  into 
usage.  Further,  for  k=4  ,  expression  (47)  holds  for  x.^  and  x^^  , 

i.e. , 


(180-0)  >  (200-30),  and 

(180-0)  >  (200-30),  respectively. 

As  shown  by  matrix  CX(I,J)  in  the  printout,  these  two  variables  are 
excluded  (fixed  at  0)  by  the  capacity  rule.  Since  the  capacity  rule 
results  in  fixing  at  least  one  variable  in  the  first  cycle,  another 
cycle  is  repeated  as  displayed  in  the  printout.  The  second  cycle 
does  not  fix  any  more  variables.  Vector  STX  is  updated  accordingly. 

Step  3:  Solve  the  relaxed  problem  (PI^)  . 

6^.7  represented  by  FLB(K)  ,  c^^  represented  by  matrix  C(I,J), 
and  solution  Y  represented  by  S0LX(J)  are  displayed  in  the  printout. 

Step  4:  Compute  LOWB. 

L0WB,  from  the  printout,  is  equal  to  749011  4375. 

Step  5:  Compare  LOWB  with  BUBS. 

Since  LOWB  <  BUBS,  go  to  Step  6. 

Step  6:  Check  if  solution  X  is  feasible  in  (P) . 

In  the  printout,  for  k=4  ,  NSUMD  =  290  >  200  ,  i.e.,  expression 
(42)  is  not  satisfied,  and  we  go  to  Step  10. 

Step  10:  As  shown  by  vector  FIX(J),  left  branching  is  possible 
and  we  go  to  Step  11. 
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Step  11:  Apply  the  bounding  rule  and  the  branching  rule. 

As  displayed  by  matrix  CX(I,J)  in  the  printout,  the  bounding 
rule  results  in  fixing  X2^  and  x^^  at  0.  Now,  for  j=4  ,  except 

for  x^  ,  all  the  x„  variables  are  fixed  at  0;  therefore  , 

is  fixed  at  1.  This  is  reflected  by  matrix  CX(I.J),  and  vector  FIX(J). 
Vector  STX  is  updated  accordingly. 

The  branching  rule  selects  x  ^  as  the  next  branching  variable. 

This  is  shown  in  the  printout  by  BR1,  and  vector  STX  is  updated 
accordingly. 

We  now  move  to  Node  3. 


Node  3: 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
printout. 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l,2, ... ,8  . 


The  facility  usage  rule  forces  facilities  1  to  5,  and  8  into 
usage.  For  k=4  ,  the  capacity  rule  excludes  x^,  and  x,^  >  i.e., 

fixes  them  at  0;  and  for  j=3  ,  all  but  being  fixed  at  0,  X23 

is  fixed  at  1.  This  is  displayed  in  the  printout  by  matrix  CX(I,J) 
and  vector  FIX(J).  Vector  STX  is  updated  accordingly. 


Although  the  capacity  rule  has  fixed  at  least  one 


ij 


variable 


during  the  initial  cycle,  another  cycle  is  not  necessary, 
by  vector  FIX(J)  which  represents  set  W,  since  we  have  an 


as  displayed 


ij 


variable 


fixed  at  1  for  each  of  the  n  columns  (activities) . 


Step  3:  Solve  the  relaxed  problem  (PR^)  • 

SOLX(J)  displays  the  solution  for  the  relaxed  problem. 


Step  A:  Compute  LOWB. 


LOWB,  shown  in  the  printout,  is  equal  to  779502.0  . 
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Step  5:  Compare  LOWB  with  BUBS. 

Since  LOWB  =  BUBS,  go  to  Step  12. 

Step  12:  Backtrack. 

We  backtrack  by  moving  to  the  parent  Node  2,  and  branching  to  the 
right  by  setting  =  0  (since  the  right  branch  has  not  yet  been  explored). 

In  the  printout,  this  is  accomplished  by  observing  the  last  entry  in  vector 
STX,  and  moving  backwards,  one  entry  at  a  time,  until  we  find  a  positive 
entry  without  1,000,000  added  to  it.  The  corresponding  x„  variable 

is  fixed  at  0,  and  we  move  to  the  right  branch  node.  Matrix  CX(I,J),  vector 
FIX(J)  and  vector  STX  are  updated  accordingly.  As  displayed  in  the  print¬ 
out,  entry  401  in  vector  STX  is  such  an  entry,  and  variable  x^  is  fixed 

at  0  for  branching  to  the  right.  This  is  shown  in  the  printout  by  BR0  as  401. 
The  updated  vector  STX  is  also  displayed. 

We  now  move  to  Node  4. 


Node  4 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  art  displayed  in  the 
printout . 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l ,2, ... ,8  . 


As  displayed  in  the  printout ,  for  k=4,  MINST=230  >  200  ,  i.e., 
the  right-hand  side  of  inequality  (47),  (s^  -  Z  dj^)  <  0  ,  and 

according  to  our  backtracking  rules,  we  backtrack,  i.e.,  go  to  Step  12. 


Step  12:  Backtrack. 


We  backtrack  to  the  parent  Node  2,  and  since  the  right-hand  branch 
has  already  beer  explored,  backtrack  to  Node  1  and  to  the  right-hand 
branch  by  fixing  x^  to  0  .  This  is  shown  in  the  printout  by  BR0  as  402, 

and  vector  STX  is  updated  accordingly. 

We  now  move  to  the  next  node,  i.e..  Node  5. 


Branch-and-Bound  Tree 


We  continue  the  branch-and-bound  procedure  from  one  node  to 
another  until  we  backtrack  to  the  root  node  and  find  that  the  right 
branch  has  already  been  explored.  The  procedure,  then,  terminates  and 
the  solution  corresponding  to  the  best  upper  bound  is  the  optimal  solu¬ 
tion. 

For  this  problem,  a  total  of  nine  nodes  are  explored  and  the 
optimal  value  equals  779502.0.  The  optimal  solution  is  x^=  X42= 

x23  =  x44  =  1  and  yl  =  y2  =  y3  =  y4  =  y5  =  y8  "  1  ‘  ThiS  ±S 
displayed  in  the  computer  printout  on  the  last  page  of  Appendix  B. 

Figure  7a  presents  the  branch-and-bound  tree  for  this  problem, 
and  shows  the  node  numbers,  the  bounds,  and  the  branching  variables. 

In  ordei  to  demonstrate  the  role  of  the  capacity  rule  and  the 
bounding  rule,  Figure  7b  displays  the  x„  variables  fixed  as  0  or  1 

by  these  rules  for  this  test  problem. 

The  cumulative  effect  of  the  various  rules,  including  the  facility 

usage  rule,  the  capacity  rule,  and  the  bounding  rule,  , lakes  the  branch- 

and-bound  procedure  quite  efficient.  Further,  the  sto  age  and  updating 

of  the  x  variables  fixed  at  0  or  1  is  done  in  a  manner  that  makes 
ij 

utmost  use  of  ti.e  relevant  information  at  the  preceding  node. 
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Figure  7a 

Branch-and -bound  tree  for  a  test  problem 
(Test  Problem  with  m=5,  n=4,  and  p=8) 


Node 

Capacity  Rule 

Bounding  Rule 

1 

x13=°  .  x33=0  ,  x14=0 

2 

X34=0  ’  X54=0 

x21=°  ,  x24=0  *  x44=1 

3 

X43-O  y  X,-  j~0  *  *23-! 

5 

X43=°  *  X34=° 

6 

XU=0  ,  x?1=0  ,  x31=0  ,  x22=0  ,  x54=0 

8 

x44  =  1 

9 

x53= ' ’  x44=0,  x54=0,  x24”' 

Figure  7b 

Variables  fixed  by  the  capacity  rule  and  the  bounding  rule 


5.  COMPUTATIONAL  TEST  RESULTS 


The  computer  program  ZIPCAP  has  been  tested  on  several  problems. 
Although  primarily  designed  for  capacitated  problems',  (i.e.,  where  the 
capacity  constraints  are  active),  the  program  can  also  be  used  for  solving 
uncapacitated  problems  as  a  special  case.  Since  the  dita  available  for 
capacitated  problems  were  limited,  some  uncapacitati-.d  problems  were  also 
considered  for  testing  the  program.  (Most  of  the  data  were  furnished  by 
Professor  Pinkur  and  are  related  to  his  work  on  multi-echelon  inventory 
systems. ) 

Table  4  presents  the  test  results  of  ZIPCAP.  In  order  to  verify 
the  optimal  solutions,  the  test  problems  were  also  solved  by  using  the 
0-1  integer  programming  code  RIP30C  [Geoff rion  and  Nelson  (1968)], 

In  the  table,  the  problem  size  shows  the  number  of  designs  (m)  , 
activities  (n).  ar.d  facilities  (p) .  This  is  equivalent  to  solving  a 
problem  having  mn+p  variables  and  n+p  constraints.  The  elapsed  time 
represents  the  time  in  seconds  to  solve  the  problem,  excluding  the  time 
to  read  and  write  the  input  and  to  write  the  outpu.:.  The  total  number 
of  nodes  explored  by  ZIPCAP  for  a  specified  set  of  options  is  also  shown. 

Both  RIP30C  and  ZIPCAP  were  run  on  an  IBM  3031  at  The  George 
Washington  University.  The  last  problem  in  the  table  was  not  run  using 
RIP30C  because  of  the  code's  capacity  limitation  to  SO  variables  and 
50  constraints. 

The  test  problem  with  m=3  ,  n=4  ,  and  p=5  has  three  variations, 
using  different  values  for  the  facility  capacities.  The  data  for  the 
variable  costs  a^j  ,  fixed  cost  b^  ,  and  the  capacity  requirements 

d...  are  given  in  the  other  document,  i.e.,  Chhabra  and  Soland  (1980). 

For  the  test  problem  with  m*5  ,  n=4  ,  and  p=8  .  runs  4a,  4b, 
and  4c  are  the  satae  except  for  nc  different  intermediate  steps'  option 


RIP30C  ZIPCAP 


19.397  277 


(ISTEP)  and  this  results  in  slight  differences  in  the  time  taken  to  solve 
the  problem.  Runs  Aa  and  4d  differ  in  that  Ad  does  nat  use  the  capacity 
rule;  the  resulting  difference  in  the  total  number  of  nodes  explored 
to  reach  the  optimal  value  points  to  the  effectiveness  of  the  capacity 
rule  in  conjun’tion  with  the  bounding  rule. 

Run  5  shows  the  results  for  an  uncapacitated  problem  with 
m=10  ,  n=8  ,  and  p=8  .  Option  ICAPR  is  not  usid  since  the  capacity 
rule  is  not  useful  for  an  uncapacitated  problem. 

Another  uncapacitated  problem  with  m=10  ,  n-30  ,  and  p=8  is 
solved  in  runs  6a  and  6b.  In  run  6a,  the  epsilon  value  (EPS)  is 
specified  as  0.002.  The  solution  value  found  by  exploring  125  nodes 
may  be  subopt  imal  but  is  guaranteed  to  be  within  -fO.2  percent  of  the 
optimal  solution  value.  Run  6b  is  made  with  an  epsilon  value  (EPS) 
of  0.0,  and  the  optimal  solution  value  is  found  in  277  nodes.  A  comparison 
of  runs  6a  and  6b  shows  that  the  number  of  nodes  ;.s  less  than  half  for  a 
solution  value  that  may  be  suboptimal  but  very  close  to  the  optimal  solu¬ 
tion  value  as  compared  to  the  number  of  nodes  for  an  optimal  solution 
value. 

In  general,  a  small  difference  between  a  solution  value  that  may 
be  suboptimal  ani  the  optimal  solution  value,  translates  into  a 
significant  difference  in  the  corresponding  number  of  nodes  and  the 
solution  time  required. 


6.  FURTHER  CONSIDERATIONS 


It  was  mentioned  in  Chapter  2  that  it  is  possible  to  consider 
alternative  formulations  of  problem  (P),  and  also  to  c  msider  choices  of 
Lagrange  multipliers  other  than  =  b^  with  the  purpose  of  obtaining 

"tighter"  bounds  which,  in  turn,  would  further  improve  the  efficiency 
of  the  branch-and-bound  procedure.  These  aspects  will  be  discussed  in 
this  Chapter. 


6.1  Alternative  Formulations 

Problem  (P)  can  be  reformulated  by  adding  additional  constraints 
such  that  the  corresponding  Lagrangian  relaxation(s) ,  if  solved,  would 
provide  "tighter"  bounds.  If  such  a  relaxation  does  not  possess  the 
integrality  property,  then  it  provides  an  equal  or  better  bound  compared 
to  that  from  an  LP  relaxation,  as  mentioned  in  Chapter  2. 

Two  alternative  formulations  of  problem  (P),  along  with  their 
Lagrangian  relaxations,  are  given  below. 
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Since  e^k  *  1  or  0,  each  constraint  of  (52)  is  either  equivalent  to 

_<  yk  (if  e^k  =  1)  or  else  is  redundant  (if  =  0)*  Problem  (API) 

thus  has,  at  nrfust,  mnp  additional  constraints  relative  to  problem  (P) . 

Two  Lagrangian  relaxations  are  now  considered  for  problem  (API). 

The  first  Lagrangian  relaxation  is  obtained  with  respect  to 
constraints  (5')  by  introducing  nonnegative  Lagrange  multipliers 
vk  ^  0  to  get 


Minimize 


E  E  a  x .  .  +  2  b  y  -  E  v  ( y,  -EE  r  ,  x..\ 
±  j  iJ  iJ  k  k'k  k  k  V  k  i  j  ijk  ) 

subject  to  (2),  (52),  and  (6),  or  equivalently. 


(ALRlv)  { 


Minimize 

E  E  x.  .  /  a  +  E  v  r  \ 

i  j  ij  ^  ij  k  k  ijkj 

-  Ek  yk  (\  - \) 

(53) 

subject  to 

Ex  =1 
i  « 

Vi 

(2) 

eikxijl  yk 

Vi, j ,k 

(52) 

x.,  ,  y  =  0  or  1 

Vi, j ,k 

(6) 

Another  Lagrangian  relaxation  of  problem  (API)  is  obtained  with 
respect  to  constraints  (5')  and  (52)  by  introducing  noni egative  Lagrange 
multipliers  v  and  A  >  respectively,  to  get 

K.  IJ  K 


Minimize  E  E  a.,  x..  +  E  b,  y, 
ij  1J  1J  k  k  k 


-  E  v. 


(yk  ■  i]  ru  k  *«) 


-  j  “  xi]k  {\  -  en>  \j) 


subject  to  (2)  and  (6),  or  equivalently. 
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(ALR1  ,)  -f 

V,  A 


Minimize 


E  Ex..  /a .  .  +  E  v  r.  +  E  e  *..,1 

ij  M1J  k  k  ijk  k  :k  ijkJ 


Subject  to 


z  yk 
k  k 

(vk  +  z  z  V 

A  1  J 

(54) 

I  X  .  . 

=  1 

vi 

(2) 

i  1J 

xij  ’ 

y^  =  0  or  1 

Vi,  j  ,k 

(6) 

For  this  problem,  the  solution  is: 


and  x  . 
1J 


y  =  0 

if  /v 

+  E  Z  X.  -  b,  \ 

1  0  * 

l  k 

i  j  1Jk  V 

=  1 

lf  (a 

+  E  E  X.  -  b,  ) 

i  .  ijk  k) 

1  V 

o 

if  i  minimizes 

1  £  j  k  k  £jl< 

+  k 

over  £ 


We  need  good  choices  of  Lagrange  multipliers  with  which  t;o 

solve  problem  (ALRlv) ,  and  of  Lagrange  multipliers  v^  and  X^k  with 
which  to  solve  problem  (ALR1  ,  ).  Problem  (ALR1  )  d^es  not  possess  the 

V ,  A  V 

integrality  property,  thus  offering  the  hope  of  a  tigr.t  bound,  but  has 
more  constraints  and  is  difficult  to  solve  compared  to  problem  (ALRlv 

which,  on  the  other  hand,  involves  more  Lagrange  multipliers. 


6.1.2  Alternative  Formulation  2 

Another  formulation  of  problem  (P)  is  similar  to  problem  (API) 
except  for  a  modification  in  constraints  (5'),  i.e.. 
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Minimize 


Subject  to 


(AP2)  M 


V 


aij  xij  +  l  Vk 

(4) 

X.  .  =  1 

IJ 

Vj 

(2) 

r  .  ..  x,  .  <  1 

ij  k  ij  - 

Vk 

(55) 

e x .  .  <  y, 

lk  ij  -  J  k 

Vi, j  ,U 

(52) 

x.  .  ,  y,  =  0  or  1 
ij  Jk 

Vi,  j  ,k 

(6) 

A  Lagrangian  relaxation  with  respect  to  constraints  (55)  and  (52)  is 
obtained  by  introducing  nonnegative  Lagrange  multipliers 
to  get 


and  X 


ijk 


Minimize  Z  I  a,,  x..  +  Z  b,  y, 
ij  i]  1J  k  k  k 


-  £  v, 

k 


\  ”  riik  \i) 

Subject  to  (2)  and  (6),  or  equivalently, 


(ALR2^)* 


Minimize 


Jk) 


Z  Z  x..  (a,,  +  Z  v  r  +  I  e  A 

t  j  iJ  \  iJ  k  k  ijk  k  lk  ijk 


+  V'  - 1  5  V) 


Subject  to 


Lx..  =1 


x  ,  y  =  0  or  1 


1  Vk 
k  K 

YJ 

Vi,  j  ,k 


(56) 

(2) 

(6) 


For  this  problem,  the  solution  is: 


yk  =  0  if  Z  l  Xijk  <_  bk  , 


i  j 

z  >: 

i  j 


*  Xijk  1  bk  * 


1  if 
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and  x, .  =  1  If  i  minimizes 
ij 


(*y  +  £  vk  ryk  +  '  e«k  x«jk  ) over  i 


Here  again,  we  need  good  choices  of  the  Lagrange  multipliers 

and  X...  with  which  to  solve  problem  (ALR2 

ijk  v  v  ,X 

6.1.3  Choice  of  Lagrange  Multipliers 

Each  of  the  relaxations  (ALRl  . )  and  (ALR2  .)  involves  p  v, 

v,A  v,A  k 

Lagrange  multipliers  and  ranp  X . . ,  multipliers.  If  we  have  good  choices 

ij  k 

of  these  multipliers,  the  resulting  solutions  of  the  relaxed  problems 

should  provide  "tighter"  bounds  (because  of  the  additional  constraints) 

than  the  bound  from  relaxation  (LR  ).  Since  relaxa;ions  (ALRl  -,  )  and 

v  v,X' 

(ALR2  ,)  are  similar  to  a  great  extent,  only  the  relaxation  (ALRl,  ■.) 

V,  A  V) A 

will  be  considered  for  further  discussion. 

By  looking  at  expression  (54)  of  the  formulation  (ALRl  ,),  a 

v,  A 

meaningful  choice  of  the  Lagrange  multipliers  and  ^  appears 

to  follow  from  setting 


v.  +  I  Z  X.  =  b,  Vk 

k  i  j  ijk  k 

so  that  each  of  the  X,.,  can  be  chosen  as 

ijk 


bk  "  Vk 
lijk  n/E  e.,\ 


if  eik  =  1  . 


0  otherwise 


The  solution  for  problem  (ALRl  .)  is  then  to  select,  from  each  column  j 

v,  A 


an  x  variable  which  minimizes 


(3ij  +  k  Vk  rijk  +  £  Hk  Vjk) 


over  i 


Arbitrary  values  were  considered  for  the  v^  (e.g.,  equal  to 

3/4  b^  ,  1/2  b^  ,  1/4  b^  .  and  0)  ,  the  X^^  were  then  computed  from 


l  ■  -*  J.V  *. 
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(58),  and  the  test  problem  with  three  designs  (m) ,  four  activities  (n) 
and  five  facilities  (p)  was  solved.  Three  cases  with  different  capacities 
s^  (as  specified  in  Chapter  5,  Table  A)  were  tried  for  the  solutions  at 

the  initial  node.  The  results,  however,  were  not  conclusive  in  terms  of 
providing  a  meaningful  choice  of  the  Lagrange  multipliers  v^  (and  of 


Since  the  relaxation  (ALR1  ,)  possesses  the  integrality  property, 

v,A 

a  choice  of  the  multipliers  as  the  optimal  values  of  the  dual  variables 
of  the  corresponding  linear  programming  problem  would  provide  a  solution 
as  good  as  the  LP  solution  (as  stated  in  Chapter  2) .  Ue  do  not  propose  to 
solve  linear  programs  as  a  part  of  our  branch-and -bound  methodology. 
However,  we  have  made  some  LP  runs,  basically  to  see  if  the  results  provide 
insight  leading  co  the  choice  of  the  Lagrange  multipliers,  and  also  to 
see  if  the  resulting  LP  solutions  are  "close"  to  the  integer  solutions. 
These  results  are  given  below. 


The  LP  formulation  (API)  corresponding  to  problem  (API)  is: 


Minimize 


Subject  to 


V 

The  constraints  x.. 

ij 


l  £ 

i  j 

I 

i 

I  £ 


a.  .  x. .  - 
ij  il 

"  I  Vk 

k 

(A) 

x.  .  =  1 

ij 

Vj 

(2) 

r  . x  ,  . 
ijk  ij 

iyk 

Vk 

(5') 

eik  xij 

<  yk 

Vij.k 

(52) 

yk 

<  1 

Vk 

(15) 

*ij  ,  yk 

>  0 

Vi, j  ,k 

(16) 

are  implicit  in  constraints  (2). 


p=5 


Problem  (API)  was  solved 
and  three  different  cases 


for  the  test  problem  with 
for  the  capacities  s^  (as 


m=3,  n=4, 

specified 


and 

in 
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Chapter  5,  Table  4).  Each  case  was  solved  using  the  IMFL  (International 
Mathematical  and  Statistical  Library)  Code  ZX3LP  on  an  IBM  3031  at  The 
George  Washington  University. 

Note  that  the  formulation  (API)  has  up  to  mnp  more  constraints  than 
the  LP  formulation  (P)  given  in  Chapter  2.  For  our  test  problem,  this 
translates  into  solving  a  problem  of  17  variables  and  50  constraints 
corresponding  to  formulation  (API)  as  against  17  variables  and  14  con¬ 
straints  corresponding  to  formulation  (P). 

Table  5  lists  the  solution  values  for  each  of  the  three  cases  with 
different  capacities  for  the  small  problem  with  throe  designs,  four  activi¬ 
ties  and  five  facilities.  The  solutions  to  problem:  (P)  and  (API), 
obtained  from  ZX3LP,  show  the  optimal  solution  values,  the  optimal  values 
of  the  variables  x„  and  y^  ,  and  the  optimal  values  of  the  dual 

variables  corresponding  to  the  Lagrangian  relaxations  (LR  )  and  (ALR1  •> ) , 

V  V  ,  A 

i.e.,  v^  associated  with  the  capacity  constraints  (j1)  and  A^.^ 

associated  with  the  constraints  (52).  The  table  also  shows  Z(P),  and 
the  Lagrangian  solution  value  Z(LR^)  obtained  by  setting  v^  =  b^  for 

all  k  at  the  root  node,  i.e.,  Z(LR,  ) . 

b 


As  expected,  the  LP  solutions  for  each  of  the  three  cases  show 
Z(AP1)  to  be  considerably  higher  than  Z(P),  and  closer  t~  Z(P),  thereby 
providing  a  tighter  bound.  As  for  the  Lagrange  multiplieis  v^  and 


A 


ijk 


the  following  relationships  are  observed. 


Z  E 

i  j 

A.  -V 

1  bk 

V k  ,  and 

+ 

> 

Z  Z 

i  j 

^ijk  —  bk 

Vk  . 

Also , 

for 

< 

II 

0  ’  Z  Z  Aijk  = 
i  J  J 

b,  ,  and 
k 

for 

b*  •  1 5  ' 

J* 

> 

o 
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Although  the  relationship  among  various  A...  values  is  not  apparent, 

lj  k 

the  above  observations  are  useful  in  further  exploring  some  good  choices 

of  the  Lagrange  multipliers  for  the  relaxation  (AI.R1  ),  as  discussed 

v ,  A 

earlier. 

As  for  the  "closeness"  of  the  LP  solution  to  that  of  the  integer 

solution,  most  of  the  solution  values  x  and  y  of  problem  (API) 

ij  k 

are  fractional,  and  their  rounding  off  to  0  or  1  does  not,  in  general, 
seem  to  correspond  to  the  optimal  integer  solution  values  x^ .  and 

yk  of  problem  (P) . 

Table  5  also  displays  ZCLR^)  at  the  root  node  for  each  of  the 
three  cases.  For  s^  =  3000  ,  Z(LR^)  =  Z(P)  ,  and  the  Lagrange  multi¬ 
pliers,  as  reflected  by  the  values  of  the  dual  variables  of  problem  (P) , 
are  equal  to  b^  for  all  k  .  This  is  expected  from  Theorem  3  and  our 

discussion  of  the  integrality  property  in  Chapter  2.  Further,  the  dual 
variables  of  (P)  for  the  first  two  cases  (i.e.,  =  7C0  Vk  ,  and 

s  =  400,...)  have  values  v,  >  b,  from  Theorem  1. 

The  Z(LR^)  values  in  Table  5,  however,  take  no  consideration  of 
the  capacity  rule  and/or  the  facility  usage  rule  of  our  branch-and- 
bound  procedure.  These  rules,  by  fixing  certain  x_  values  to  0 

or  1,  and  by  forcing  certain  facilities  into  the  solution,  provide 

an  improved  lower  bound.  As  per  our  branch-and-bound  procedure  the 

improved  lower  bound  at  the  root  node  is  obtained  by  solving  problem 

(PR,).  For  example,  for  s.  =  700  Vk,  the  values  of  Z(LR,)  and  Z(PR. ) 

-L  k  b  1 

are  shown  in  Figure  8.  The  figure  also  shows  the  values  of  Z(P) , 

Z(APl),  and  Z(P).  The  branch-and-bound  procedure  rules  provide  an 

improved  value  of  the  lower  bound  Z(PR^)  compared  to  Z(P).  It  appears 

that  some  good  values  of  the  Lagrange  multipliers  of  the  relaxation 

(ALRl  ,),  if  found,  could,  in  conjunction  with  these  rules,  provide 
v ,  A 

significant  improvement  over  Z(AP1),  and  without  the  need  to  solve  an 
LP  problem. 


" .. ii w woawt 
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6.2  Subgradient  Method 

It  was  mentioned  in  Chapter  2  that  setting  the  Lagrange  multipliers 
equal  to  for  all  k  provides  a  good  starting  point  in  solving 

the  Lagrangian  relaxation  (LRv)  of  problem  (P) .  From  Theorem  3,  this 

choice  is  optima]  (in  terms  of  providing  the  tightest  lower  bound)  if 
the  resulting  solution  is  feasible  in  problem  (P) .  In  other  cases, 
i.e.,  where  the  resulting  solution  is  not  feasible  in  problem  (P) ,  this 
choice  is  generally  not  optimal  and  it  is  possible  to  tighten  the  bounds 
by  considering  values  of  v^  b^  (from  Theorem  1).  Oi.e  method  that 

seems  useful  in  providing  such  a  choice  is  the  subgradient  method. 

The  subgradient  method  is  an  adaptation  of  the  gradient  method 
in  which  gradients  are  replaced  by  subgradients.  Through  a  heuristic 
choice  of  the  step-size,  this  method  has  been  successfully  used  to 
provide  improved  bounds  and  sometimes  optimal  solutions  [for  details 
see  Held,  Wolfe,  ar.d  Crowder  (1974),  Fisher  (1978),  and  Christofides 
(1980)].  The  fundamental  theoretical  result  is  that 

Z(LR  8)  Z(D)  if  t8— *0  and  Z  t9— ►  «>  as  g  — *•  , 

V  q=0 

where  t8  is  the  positive  step-size  t  for  the  g'.h  iteration, 
and  Z(LR^8)  is  the  solution  value  of  the  relaxed  problem  (LRv)  using 

v^  values  obtained  at  the  gth  iteration. 

e+1 

In  the  case  of  problem  (P) ,  the  step-size  t6  ) or  iteration 
g  +  1,  given  that  we  have  a  solution  of  (LRv8),  is  given  by 


t-8+1 


A8+1[Z*  -  Z(LR  8)] 


'  l  j  rijk  Xij 


8  I  I  2 


(59) 


g+1 

where  A  is  a  scalar  and  generally  between  0  and  2,  and 

Z  is  an  upper  bound  on  Z(LRv8),  frequently  obtained  by  applying 


a  heuristic  to  solve  problem  (P) . 
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Given  the  vector  of  multipliers 


.B+l 


is  generated  by 


z  z  r 

i  J 


ijk 


(60) 


where  we  enforce 


in  our  case  of  problem  (P)  (because  of 


Theorem  1) . 


Justification  for  these  rules  and  computational  results  of 
applications  of  the  subgradient  method  are  given  in  Held  e_t  al^  (1974). 

The  scalar  A  is  generally  initiated  by  setting  A^  =  2  and  halving 
subsequent  A's  whenever  the  resulting  solution  value  has  failed  to 
increase  in  some  fixed  number  of  iterations.  This  rule  has  performed 
well  empirically  [Held  et  al  (1974)  and  Fisher  (1978)]. 

For  the  test  problem  with  three  designs,  four  activities,  five 
facilities,  and  the  capacities  s^  =  400,  S2  *  400,  =  1000, 

=  400,  s^  =  400,  the  Lagrangian  solution  obtained  it  the  root  node 

by  setting  =  b^  for  all  k  ,  i.e,,  the  solution  to  problem  (PR^) 

is  infeasible  ii  problem  (P) ,  i.e.,  it  violates  the  capacity  constraints. 
It  seems  that  the  subgradient  method  could  be  uselul  in  considering 
v  ^  b  with  the  ultimate  purpose  of  obtaining  a  tighter  lower  bound, 

and,  depending  on  the  revised  sol ution(s) ,  possible  jmprovement  in  the 
best  upper  bound.  Another  possiblity  could  be  to  first  arbitrarily 
increase  the  relevant  v^  values  by  a  small  percentage  of  the  b^ 

values  and  then  solve  problem  (LR^),  hopefully  to  improve  the  lower 

bound;  and  thereafter  to  use  the  subgradient  method  for  obtaining 
subsequent  values  of  v^  ,  and  tightening  the  bounds. 


Both  of  the  areas  discussed  above,  i.e.,  the  consideration  of 
alternative  formulations  of  problem  (P),  and  the  application  of  the 
subgrad iont  method,  and  their  combination,  seem  useful  for  continued 
research  in  terms  of  further  improving  the  branch-ani-bound  procedure 
for  solving  the  multiactivity  multifacility  capacity-constrained  0-1 
assignment  problems. 
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FORTRAN  IV  G  LEVEL  21 


MAIN 


DATE  =  80315 


11/07/24 


0001 


0002 

0003 

COOL 

0005 


0006 
000  7 


0008 

0009 

0010 

0011 

0012 

0013 


00  lL 


0015 
0016 
00  17 
0018 
0019 
0020 
0021 
0022 


C 

C 

C 

C 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 


c 

c 

c 


c 


c 


c 

c 

c 

c 


ZIPCAP,  A  ZERO-ONE  INTEGER  PROGRAM  IS  DESIGNEO 
TO  SOLVE  MULTIACTIVITY  MULT  I F AC  I L I TY  CAPACITY- 
CONSTRAINED  PROBLEMS  HAVING  VARIABLE  AND  FIXED 
COSTS.  IT  ALSO  SOLVES  UNC Al  AC  I TATEO  PROBLEMS  AS  A 
SPECIAL  CASE 

INTEGER  0(35,35,301,  A<35,35),  CX(35,3f),  E135.30), 

1  8(301,  8S0LX (35),  BSQLY(30),  FL8 ( 30 1 , F I  X ( 35 ) , F I  X I ( 3 5  I  , 

2  FUB ( 30  1 ,  S ( 30 1  ,  S0LX(35),  STX(122r» 

REAL  M I NC ( 3  5 1 ,  NMINC<35) 

DIMENSION  0(35, 35),  0IFBRO5),  KT2I351,  MIND(3>) 

-INTEGER  BRO,  BR1,  FC,  FCUB,  P 

REAL  LOWB,  MAXDIF,  MINSC 

»*»***» **»****QPT IONS  AVAILABLE:  IINPT,  ICAPR,  1STEP,  IUNCAP.EPS 
1 INP T  =  1  IF  INPUT  LISTING  DESIRED;  0  OTHERWISE 
ICAPR  =1  IF  CAPACITY  RULE  TO  BE  USED;  0  OTHERWISE 
ISTEP=0  IF  LISTING  OF  INTERMEDIATE  STEPS 
NOT  DESIRED.  ISTEPM  IF  SUMHARY  OF  BRANCH  L 
BOUND  NODES  DESIRED.  ISTEP=2  IF  DETAILED 
LISTING  Of  INTERMEDIATE  STE^S  DESIRED. 

IUNCAPM  IF  SOLVING  AN  UNC  Ai' AC  I T  ATED  PROBLEM, 

0  OTHERWISE. 

EPS=  A  FRACTIONAL  VALUE  IF  SUBOPTIMAL 
SOLUTION  DESIRED,  E.G.,  EPSILON  AS  0.005 
IMPLIES  SOLUTION  TO  BE  WITHIN  vC.5  PERCENT 
OF  THE  OPTIMAL  SOLUTION.  EPS=0.0  IF  OPTIMAL 
SOLUTION  DESIRED. 

ET  =  ELAPSED  TIME  IN  SECONDS,  IF  SPECIFIED,  AT 
WHICH  THE  NOOE  AND  BOUNDS  RELATED  INFORMATION 
IS  PRINTED.  THIS  IS  USEFUL  IN  \  SITUATION  IF 
I STEP=0  AND  THE  PROGRAM  TERMINATES  before 
REACHING  THE  FINAL  SOLUTION. 

•  ***>  ««»»»»»»» RE  AD  INPUT  DATA*************-****«*************** 

READ  10,  IINPT,  ICAPR,  ISTEP,  IUNCAP,  EPS,  ET 
10  FORMAT  (Ml,  F6.5,  F10.3) 

M*  NUMBER  OF  DESIGNS 
N=  NUMBER  OF  ACTIVITIES 
P=  NUMBER  OF  FACILITIES 

READ  20 ,M, N, P 
20  FORMAT  (3151 

A ( 1 , J  )  :  VARIABLE  COST  MATRIX 
READ  30,  ((AII.J),  1=1,M),J=1,N) 

30  FORMAT  (8110) 

B(K) :  FIXED  COST  VECTOR 
READ  30,  <B(K) , K  =  1 , P ) 

IF  (IUNCAP. EQ.l)  GO  TO  40 

S ( K I :  CAPACITY  LIMIT  VECTOR;  REQUIRED  ONLY 
IF  1 UNC  A  P  =  0 
READ  -.0,  (SIK),K=1,P) 

0(1, J,K):  CAPACITY  USAGE  MATRIX;  REQUIRED 
ONLY  IF  IUNCAP=0 

DO  32  K= 1 , P 

READ  30 , ( ( D(  I, J,K I  ,  IM,Ml,J=l,NI 
32  CONTINUE 
DO  37  KM.P 
00  37  1  =  1, M 

IF  (0(1,1, Kl.EQ. 01  GO  TO  35 
£  (  I, KIM 
GO  TO  37 


00000010 

00000015 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

00000080 

00000090 

000001)0 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000253 

00000256 

00000260 

00000263 

00000266 

00000270 

00000273 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00000340 

00000350 

00000360 

00000370 

00000380 

00000390 

00000400 

00000410 

00000420 

00000430 

0000C440 

00000450 

00000460 

00000470 

00000430 

0000049C 

00000500 

00000510 

00000520 

00000530 
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0023 

35 

E (  I ,  K I *0 

00000540 

0024 

37 

CONTINUE 

00000550 

i 

0025 

GO  TO  90 

00000560 

C 

E ( I ,K  1 :  DESIGN-FACILITY  MATRIX!  REQUIRED  ONLY 

00000570 

C 

IF  I UNCA  P=1 

00000580 

C 

0026 

40 

RE  AO  45,((E(I,K),l*l,M),K=l,P> 

00000590 

0027 

45 

FORMAT  (8011) 

00000600 

0028 

DO  30  K=1,P 

00000610 

o 

0029 

5  ( K  )  =N 

00000620 

0030 

DO  75  1=1, M 

00000630 

0031 

IF  (Ed.KI.EQ.l)  GO  TO  65 

00000640 

© 

0032 

00  60  J=1,N 

00000650 

0033 

01 I, J,K)=0 

00000660 

0034 

60 

CONTINUE 

00000670 

o 

0035 

GO  TO  75 

00000680 

0036 

65 

DO  70  J*l,N 

00000690 

0037 

D(  I , J ,K ) =1 

00000700 

o 

0038 

70 

CONTINUE 

00000710 

0039 

75 

CONTINUE 

00000720 

0040 

80 

CONTINUE 

00000730 

0 

C 

00000740 

0041 

90 

PRINT  95,  IINPT,  ICAPR,  ISTEP,  IUNCAP,  EPS,  ET 

00000750 

0042 

95 

FORMAT  (  •!•,•  OPTIONS  SELECTED  :  IINPT*', 11, 

00000760 

0 

1 

L  •  I C  APR* ‘,11,  •  ISTEP**, 11,  •  1  UNCAP  =  • ,11, 

00000770 

2  •  EPS=‘,F8.5,  •  ET=‘ ,  F10.3///1 

00000780 

0043 

IF  (IINPT. EQ.O)  GO  TO  168 

00000790 

c 

0044 

PRINT  100,M,N,P 

00000800 

0045 

100 

FORMAT  (‘O',  T55 ,  ‘INPUT  OATA‘,/lX,  T55,  //////IX , 

00000810 

1T41,  ‘NUMBER  OF  DESIGNS  (M)*»,  *X ,14//1X,T*1 , 

2  *  NUMBER  OF  ACTIVITIES  (N)**,  IX.I4//1X,  T<.  1 , 

3 • NUMBER  OF  FACILITIES  <P)=*,IX,  I4///1 
PRINT  105 

105  FORMAT  (  4X.  ‘VARIABLE  COST  MATRIX  A I  I , J ) • ,/4X, 

1 '  —  --  — .  — — '  ,/ ) 

DO  110  1=1, M 

110  PRINT  115,  I,  (A(I,J),J=1,N) 

115  FORMAT  l‘0‘,  T6,  'I*',  13,  4X.8I13,  4(/,  14r, 8113)1 
PRINT  120 

120  FORMAT ( *0* ,//4X, ‘FIXED  COST  VECTOR  B(K)‘,/4X, 

1  ‘ 1  ■  ■■■■■■  ■  - »,  /) 

PRINT  122,  <B(K),K=1,PI 

122  FORMAT  «'0‘,  T15,  8113,  3C/,  14X,8I13») 

PRINT  125 

125  FORMAT) ‘0‘ ,//4X,  ‘CAPACITY  LIMIT  VECTOR  S(M*,/4X, 

1»  — ------ ■—  . . . . — »,  /) 

PRINT  128,  (S(KI,K=1,P) 

128  FORMAT  l‘0‘,  T15,  3113,  31/,  14X.8U3)) 

PRINT  130 

130  FORMATC *0‘ ,//4X,  ‘CAPACITY  USAGE  MATRIX  D1I*J,K*‘,/4X, 

1‘- - - •  ,/) 

00  150  K=1.P 

PRINT  135, K 

135  FORMAT  ( ‘0 • ,//5X , ‘K  =  • , 1 3 / ) 

00  145  1=1, M 

PRINT  140, 1  ,(0( I,J,KI ,  J  =  l ,N  • 

140  FORMAT  (*0‘,  T6,  •!=',  13,  4X.8I13,  4(/,  14X.B113)) 

145  CONTINUE 

150  CONTINUE 


00000820 
00000830 
000008*0 
00000850 
00000860 
00000870 
00000880 
00000890 
00000900 
00000910 
00000920 
00000930 
00000940 
00000950 
00000960 
00000970 
00000980 
00000990 
00001000 
00001010 
00001020 
00001030 
0000 1G40 
00001050 
00001060 
OCOO 10  70 
00001080 
00001090 
00001100 
00001110 
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0069 

PRINT  155 

00001120 

0070 

155 

FORMAT! *0* ,//4X , 'DES IGN-FACIL 1TY  MATRIX  EII,K!*,/4X, 

00001130 

( 

1 ' - - — . . . . ',/) 

00001 140 

0071 

00  160  1*1, M 

00001150 

0072 

PRINT  158,  I,  (EII,K)»K*1,P) 

00001160 

( 

0073 

158 

FORMAT  ('O',  T6 ,  •!=',  13,  4X,8I13,  3«/,  14X.8I13II 

00001170 

0074 

leO 

CON .  INUE 

00001 180 

0075 

i6e 

IF  <ISTEP.EQ.01  GO  TO  190 

00001190 

3 

0076 

IF  USTEP.EQ.il  GO  TO  175 

00001200 

0077 

PRINT  170 

00001210 

0078 

170 

FORMAT  1 'O', ///55X, 'DETAILED  LISTING  OF  STEPS',/) 

00001220 

,  X 

0079 

GO  TO  190 

00001230 

0080 

175 

PRINT  180 

00001240 

0081 

180 

FORMAT  ( 'O', ///55X, 'SUMMARY  OF  STEPS',/! 

00001250 

o 

C 

00001260 

0082 

190 

BUB=9999999. 

00001270 

0083 

8U8S=  BU8/  I  1 .0  *EPS I 

00001280 

o 

0084 

NSX=0 

00001290 

0085 

NOO=1 

00001310 

0086 

IBNOO-1 

00001315 

e 

0087 

INE  T  =  0 

00001320 

0088 

IN3ET  =0 

00001330 

0089 

00  205  J  =  l,N 

00001390 

c 

0090 

f: XI J)=0 

00001400 

0091 

KT2 I J 1 “0 

00001410 

0092 

DO  205  1=1, M 

00001420 

c 

0093 

CX11,J)*0 

00001430 

0094 

205 

CONTINUE 

00001433 

0095 

LQ  1=0 

00001436 

0096 

LQ2=0 

00001440 

0097 

LR2=0 

00001443 

0098 

CALL  TIMET  1 ITO) 

000014*5 

o 

0099 

IF  USTEP.EQ.O)  GO  TO  208 

00001448 

0100 

PRINT  220, NOD 

00001450 

0101 

208 

IFINSX.EQ.OI  GO  TO  283 

00001453 

0 

C 

CX(I,J)  CONTAINS  FIXED  AND  FREE  XII, J)  VARIABLES. 

00001456 

C 

STXIINSI  CONTAINS  FIXED  XII, J)  VARIABLES. 

00001460 

C 

CXII,J)  AND  STXIINSI  ARE  UPDATFD  BY  THE  CAPACITY 

00001480 

o 

c 

RULE,  THE  BOUNDING  RULE,  AND  THE  RULE  FOR 

00001490 

c 

BRANCHING  AND  BACKTRACKING. 

00001500 

c 

IN  CXII.JI  A  FIXED  VARIABLE  <S  RECORDED  AS  1  OR 

00001505 

0 

c 

2,  AND  A  FREE  VARIABLE  AS  0. 

00001510 

c 

A  VALUE  OF  l  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 

00001515 

c 

IS  FIXED,  AND  FIXIJ)  IS  SET  EQL'A L  TO  1  IMPLYING 

00001520 

o 

c 

THAT  COLUMN  J  HAS  A  FlXEO  VARIABLE  OF  VALUE  1. 

00001525 

c 

A  VALUE  OF  2  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 

00001530 

c 

SHOULD  NOT  BE  CONSIDERED  FOR  CURRENT  COMPUTATIONS. 

00001535 

9 

c 

AN  XI  I, J)  RECORDED  IN  CXII,J)  AS  1  DUE  TO  THE 

000015*0 

c 

BRANCHING  RULE  IS  RECORDED  IN  STXI.’NS)  AS  X*100*J. 

00001545 

c 

AN  XII, J!  RECORDED  IN  CXII,J»  AS  l  DUE  TO  THE 

00001550 

(. 

c 

CAPACITY  RULE  OR  THE  BOUNDING  Rl,LE  IS  RECORDED  IN 

00001555 

c 

STXIINSI  AS  IX*100*J)*100C000. 

00001560 

c 

AN  XII, J!  RECORDED  IN  CXII.J)  AS  2  IS  RECORDED  IN 

00001565 

c 

STXIINSI  AS  -IX*100*J 1-1000000. 

00001570 

0102 

210 

IF  1  1STEP.IQ.0I  GO  TO  225 

00001580 

0103 

215 

PRINT  220, NOD 

00001590 

0104 

220 

FORMAT  1 'O', //6X, 'NODE  NUMBER',  17/) 

00001600 

C 

00001610 
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c 

BRO  IS  THE  RIGHT  BRANCHING  VARIABLE 

00001615 

0105 

225 

LXbBRO 

0000)620 

0106 

I X*LX/ 100 

00001630 

0107 

JX=LX-IX*100 

00001640 

0108 

CXI IX,JX)b2 

00001650 

0109 

KT2I JX)*KT2( JX)*1 

00001660 

0110 

FIX! JX)=0 

00001720 

0111 

LQ1bLQ1-1 

00001725 

0112 

IF  IKT2IJXI.CT.IM— 11)  GO  TO  270 

00001730 

0113 

DO  255  1=1, M 

00001740 

0114 

IF  (CXI  I,JX).EQ.2)  GO  TO  255 

00001750 

0115 

CXI  I  ,JXI«l 

00001760 

0116 

NSX=NSX»1 

00001763 

0117 

STX(NSX)=  I I*100« JXI 41000000 

00001766 

0118 

FIX! JX)»l 

00001770 

0119 

LQ l=LQl+l 

00001780 

0120 

FIX! (JX)=1 

00001790 

0121 

GO  TO  270 

0C001800 

0122 

255 

CONTINUE 

00001810 

0123 

270 

LQ2=0 

00001820 

0124 

LR2=0 

00001825 

0125 

GO  TO  283 

00001830 

0126 

272 

IF  (I  STEP. EG. 0)  GO  TO  276 

00001840 

0127 

PRINT  220, NOD 

00001850 

c 

00001853 

c 

BR1  IS  THE  LEFT  BRANCHING  VARIABLE 

00001856 

0128 

276 

LQ2«0 

00001860 

0129 

LR2=0 

00001866 

0130 

LX*BR1 

00001870 

0131 

1X-LX/100 

00001875 

0132 

JX=LX-IX*100 

00001880 

0133 

CX(IX,JX)=1 

00001885 

0134 

FIX! JX)=1 

00001890 

0135 

LQ1=L01*1 

00001892 

0136 

00  279  1  =  1, M 

00001895 

0137 

IF  UX.EQ.I)  GO  TO  281 

00001897 

0138 

279 

CONTINUE 

00001900 

0139 

281 

FIXII JXI=IX 

00001902 

0140 

283 

IF  1 1  STEP .NE .2 )  GO  TO  303 

00001905 

0141 

285 

DO  295  1=1 ,M 

00001910 

0142 

PRINT  290.  t.(CX(I,J),J=l,N) 

00001920 

0143 

290 

FORMA  T  (/5X,'CX(I,J),.4X .*!  =  '. 13. 

2  X ,  20 14/23X .  2014) 

00001930 

0144 

295 

CONTINUE 

00001940 

0145 

PRINT  297, (FIX! J),J=I,N) 

00001950 

0146 

297 

FORMAT  l/5X,*FIX(J)ftl2X,  20I4/23X, 

2014) 

00001960 

c 

»,**«*, **.*,..AppLY  CAPACITY  RULE*** 

00001970 

c 

AND  UPDATE  CXI  I.J)  AND 

ST  X 1  INS). 

00001980 

0147 

303 

DO  307  K=1,P 

00002000 

0148 

FLB I K ) =0 

00002015 

0149 

307 

CONTINUE 

00002025 

0150 

310 

DO  2000  K=1,P 

00002030 

c 

FIND  THE  SUM  OF  MINIMUM  DII.J.K)  OVER 

EACH  J  FOR  A 

00002040 

c 

GIVEN  K,  I.E.,  MINSDb 

SUM  OF  v I ND 1 J  ) 

00002050 

0151 

M I NS0=0 

00002060 

0152 

00  400  J=1,N 

00002070 

0153 

IFIFIXl Jl.EO.OI  GO  TO  350 

00002030 

c 

IF  FIX(J)=l,  SET  MtN0(J)=0(I,J,KI  FOR 

CXI !,JI=t 

00002090 

c 

ANO  MOVE  TO  NEXT  COLUMN  J 

00002 100 
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0154 

1N0I*FIXI(J1 

00002110 

0155 

M1ND(J1*0(IND1,J,K! 

00002120 

< 

0156 

GO 

TO  800 

00002130 

0157 

350 

LK=0 

00002160 

0158 

1*1 

00002170 

(. 

0159 

HINDI J)*D(I,J,K) 

00002180 

C 

SKIP  0(1. J.K)  WHEN  CX  ( I , J 1*2  C  HOVE  TO  NEXT  ROW  I 

00002190 

OloO 

400 

IF ( CX ( I  .3 1  .EO .2  I  GO  TO  600 

00002200 

0161 

500 

TF(0(I«J,K l.LT.MIND ( J 1 1  HINDI J 1*0 ( I,J,X1 

00002210 

01o2 

GO  TO  700 

00002220 

0163 

600 

LK*LK*1 

00002230 

o 

0164 

IFII.GT.LKI  GO  TO  700 

00002240 

0165 

i*i*i 

00002250 

0  166 

HINDI J 1 *01 I ,  J  ,K  1 

00002260 

o 

0167 

GO  TO  750 

00002270 

OloS 

700 

i«i*i 

00002280 

0169 

750 

IFCI.LE.HI  GO  TO  400 

00002290 

c 

0170 

800 

HI NSD=M1NS0*M1ND 1 J 1 

00002300 

0171 

900 

CONTINUE 

00002310 

0172 

910 

IF 

I ISTEP.NE.2)  GO  TO  960 

00002320 

o 

0173 

PRINT  950,  K,  M1NSD , (HINDI J) ,J*1 ,N) 

00002330 

0174 

950 

FORHAT  (•0,,,K,MINS0,(MIND(J1,J*1,NI*,10I10,4(/,44X,8I101 1 

00002340 

0175 

960 

IF  (MINSD.EC.O)  GO  TO  975 

00002342 

0 

0176 

965 

IF (  FL8 (K 1 .EQ.l 1  GO  TO  975 

00002344 

0177 

970 

F..8  (K 1  *1 

00002346 

0178 

975 

IF  (IUNCAP.EO.il  GO  TO  2000 

00002348 

o 

0179 

978 

IF  1ICAPR.EQ.0)  GO  TO  2000 

00002349 

c 

FINO  BALANCE  AVAILABLE  CAPACITY  I3ALD  FOR  A  GIVEN  K 

00002350 

c 

IF  I8AL0  IS  NEGATIVE,  THEN  BACKTRACK. 

00002360 

L 

0130 

980 

IB ALD*S IK 1— MINSO 

00002380 

0181 

1000 

IF 

(IBALD.LT. 01  GO  TO  6200 

00002390 

0182 

00 

1500  J-l.N 

00002400 

0 

c 

SKIP  COLUMN  J  IF  FIX(JI*1 

00002 410 

0183 

IF  (FlXIJI.EQ.il  GO  TO  1500 

00002420 

0184 

00 

1300  1*1, H 

00002430 

0 

c 

SKIP  ROW  I  IF  CX ( I , J 1 *2 

00002440 

0185 

1100 

IF(CX( I.J1.E0.2I  GO  TO  1300 

00002450 

c 

COHPUTE  DIFFERENCE  BETWEEN  OII.J,Kl  AND  HINDU!. 

00002470 

o 

c 

IF  IT  IS  MORE  THAN  AVAILABE  BALANCE,  SET  CX(I,JI*2 

00002480 

0186 

1200 

IDlF0=0(I,J,K)-MINDIJ) 

00002490 

0137 

IF  ( ( I DI FD— 1 BALD1.LE.0)  GO  TO  1300 

00002510 

o 

0188 

CXU,J1=2 

00002520 

0139 

NSX=NSX*1 

00002523 

0190 

STX ( NSX )=-(!* IO0* J 1-1000000 

00002526 

J 

c 

LQ2  COUNTS  THE  NUMBER  OF  CXII,J)  VALUES  SET  EQUAL 

00002530 

c 

TO  2  IN  A  CYCLE 

00002540 

0191 

LQ2*LQ2  *1 

00002550 

o 

c 

KT2IJ)  KEEPS  AN  ACCOUNT  OF  - CX ( I , J 1  VALUES  SET  EQUAL 

00002560 

c 

TO  2  FOR  COLUMN  J 

00002570 

0192 

KT2 ( J 1 =KT2 ( J 1 *1 

00002580 

1 

c 

FOR  COLUMN  J,  IF  ALL  BUT  ONE  r.X(I,Jl  VALUES  ARE 

00002590 

c 

EQUAL  TO  2,  SET  THAT  CX(!,JI*1  C  SET  F1X(J!*1 

00002600 

0193 

IF(KT2(J).LT.( M— 1 ) 1  GO  TO  1300 

00002610 

1 

0194 

00 

1250  LR=l,M 

00002620 

0195 

IF(CX(LR,J1.EQ.2I  GO  TO  1250 

00002630 

0196 

CaC  R,J  1  =  1 

00002640 

0197 

NSX=NSX*l 

00002643 

0198 

STX(NSX1=  ( LR*100*J  >  *1000000 

00002646 
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0199 

FIX( J 1*1 

00002650 

c 

LQ1  KEEPS  AN  ACCOUNT  OF  COLUMNS  FOR  WHICH  FIX(J>*1 

00002655 

0200 

LQ1=LQ1+1 

00002660 

C 

FIXI(J)  SPECIFIES  INDEX  I  FOR  WHICH  FIX(J)*1 

00002662 

0201 

FI  XI ( Jl *LR 

00002665 

0202 

GO  TO  1500 

00002670 

0203 

1250 

CONTINUE 

00002680 

020A 

1300 

CONTINUE 

00002690 

0205 

1500 

CONTINUE 

00002700 

0206 

1800 

IF  (15T£P.N£»2)  60  TO  2000 

00002710 

0207 

PRINT  1900 «  K*  LQ2>  LQ1 

00002720 

0208 

1900 

FORMAT  ( *0  * »  *K* • » 13»  *  LQ2=*,I3,  •  LQ1*1,  I3> 

00002730 

0209 

00  1930  1*1, M 

00002740 

0210 

PRINT  290,  I*(CX<I*J),J=1«N) 

00002750 

0211 

1930 

CONTINUE 

00002770 

0212 

PRINT  297,  <FIXCJI,J*1,NJ 

00002780 

0213 

2000 

CONTINUE 

00002800 

C 

A  CYCLE  EXAMINES  ALL  THE  FACILITIES. 

00002803 

C 

IF  IN  A  CYCLE,  THE  CAPACITY  RULE  RESULTS  IN  SETTING 

00002810 

C 

ADDITIONAL  CX(I,J>  VALUES  EQUAL  TO  2,  THEN  REPEAT 

00002820 

C 

THE  CYCLE.  BUT  IF  FIXIJ)*1  FOR  ALt  J,  THEN  DO  NOT 

00002830 

c 

REPEAT  THE  CYCLE. 

00002835 

0214 

IF  ILQ1.E0.N)  GO  TO  2400 

00002840 

0215 

IF  (LQ2.EQ.LR2I  GO  TO  2400 

00002845 

0216 

2200 

LR2*LQ2 

00002860 

0217 

GO  TO  310 

00002870 

c 

00002880 

c 

UPOATE  VECTOR  OF  FACILITIES  FLBCKl  FOR  COMPUTING 

00002890 

c 

C 1 X  *  J 1  MATRIX  £  COWER  80UN0.  IT  HAS  VALUE  l  IF  A 

00002900 

c 

FACILITY  IS  USED,  OTHERWISE  IT  HAS  0  VALUE. 

00002910 

0218 

2400 

DO  3000  J*1,N 

00002950 

0219 

IF  (FIxm.EQ.O)  GO  TO  3000 

00002960 

0220 

INOI*FI XI ( Jl 

00002970 

0221 

DO  2550  K*1 »P 

00002990 

0222 

IF  (EIINDI.KI.EQ.O)  GO  TO  2550 

00003000 

0223 

IF  (FLB(K).EQ.l)  GO  TO  2550 

00003010 

02  24 

FLB(KI=1 

00003020 

0225 

2550 

CONTINUE 

00003030 

0226 

3000 

CONTINUE 

00003060 

0227 

IF  ( I STEP.NE .2 )  GO  TO  3150 

00003070 

0228 

PRINT  3100,  (FLB(K) ,K=l,P> 

00003080 

0229 

3100 

FORMAT! *0*  »• (FLB (K) , K* 1 , P 1  •,  201 4/16X, 201 4 > 

00003090 

C 

COMPUTE  COST  MATRIX  C(1,J)  FOR  THE  RELAXED  PROBLEM 

00003100 

0230 

3150 

DO  3400  J*1,N 

00003110 

0231 

DO  3300  1*1, M 

00003120 

02  32 

BSUM=0.0 

00003130 

0233 

DO  3200  K*1,P 

00003140 

0234 

IF  ( FLB (K 1 .EQ .1 >  GO  TO  3200 

00003150 

02  35 

IF  (Ell  .KI.EQ.O)  GO  TO  3200 

00003160 

0236 

BS UM*BSUM+ ( 3 ( K )  •  (F10AT<0(I,J,K*)/  FL0AT1SIKII 1) 

00003170 

0237 

3200 

CONTINUE 

00003180 

0238 

3250 

C»  I,J)*A(I,J) *BSUM 

00003190 

02  39 

3300 

CONTINUE 

00003200 

0240 

3400 

CONTINUE 

00003210 

0241 

IF  ( I STEP.NE .2 )  GO  TO  3445 

00003220 

0242 

00  3430  1*1, M 

00003230 

0243 

PRINT  3420,  I,  (C< I, Jl ,J=l,N) 

00003250 

0244 

3420 

FORMAT  I/5X,  •CII,J)',5X,  •!*',  13, 2X,  5F15.4, 

00003260 

01 
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1  6I/23X,  5F15.4II 

00003265 

0245 

3430 

CONTINUE 

00003270 

C 

FIND  SUM  OF  MINIMUM  C(I,J>  VALUES  OVER 

EACH  J, 

00003290 

C 

I.E.,  MINSC  =  SUM  OF  MJNCU). 

00003300 

C 

IF  F  I  X  (  J  I  *  l ,  THEN  MINC(J)=C( J,J1  WHERE 

CX(1,J)*1 

00003310 

0246 

3445 

MINSC=0.0 

00003320 

0247 

00  3900  J*t,N 

00003340 

0248 

IF  (FIXU1.EQ.01  GO  TO  3500 

00003350 

0249 

INDI=FIXI ( J) 

00003360 

0250 

MINC( J)=Cl INDI ,J) 

00003370 

0251 

SOLX(J)=INOI 

00003380 

0252 

GO  TO  3850 

00003410 

0253 

3500 

LK=0 

00003430 

0254 

1  =  1 

00003440 

C 

SKIP  CII.JI  ELEMENT  IF  CX(!,J):2  t  MOVE 

TO  NEXT  I 

00003470 

0255 

3550 

IF  (CX< l,J J.EQ.2J  GO  TO  370: 

00003480 

0256 

IF  ((1-LKI.EQ.l)  GO  TO  3600 

00003485 

0257 

IF  (CII.JI. GE.MINCIJJ)  GO  TO  3750 

00003490 

0258 

3600 

MlNC(JI*CU,JI 

00003500 

02  59 

IMINxI 

00003510 

0260 

GO  TO  3750 

00003520 

0261 

3700 

LK=LK*1 

00003530 

0262 

3750 

1=1*1 

00003590 

0263 

3800 

IF  (I.LE.MI  GO  TO  3550 

00003600 

02o4 

SQLX 1 J  I  = I M IN 

00003610 

0265 

3850 

minsc=minsc*minciji 

00003620 

0266 

3900 

continue 

00003630 

0267 

IF  USTEP.NE.2I  GO  TO  3940 

00003640 

0260 

DO  2920  J=1,N 

00003650 

0269 

PRINT  3910,  J,MINC« J) ,SOLX(J) 

00003660 

0270 

3910 

FORMAT  (•0*,*J,MINC(JI,S0LX(JI‘,  I5,F15.4,I6> 

00003670 

0271 

3920 

CONTINUE 

00003680 

C 

COMPUTE  FIXED  COST  FC  FOR  LCWER  BOUND 

00003710 

02  72 

3940 

FC=0 

00003720 

0273 

DO  4000  K=1,P 

00003730 

02  74 

IF  (FLB(Kl.EQ.O)  GO  TO  4000 

00003740 

02  75 

3950 

FC=FC*B(K) 

00003750 

02  76 

'  000 

CONTINUE 

00003760 

c 

00003770 

0277 

4050 

L0W8=MINSC*FC 

00003780 

0278 

IF  (ISTEP.EQ.OI  GO  TO  4150 

00003790 

0279 

PRINT  4120,  MINSC,  FC,  L0W8 

00003800 

0280 

4120 

FORMAT  (•O',*  MINSC,  FC,  LOWS  F15.4,  115,  F15.4I 

00003810 

C 

COMPARE  LOWER  BOUND  WITH  BEST  UPPER  BOUNO  STAR 

00003820 

C 

BUBS  WHICH  EQUALS  BU8/tl*EPS).  IF  L0W8 

IS 

00003830 

C 

GREATER  THAN  OR  EQUAL  TO  BUBS,  THEN  BACKTRACK 

00003840 

0281 

4150 

IF  ( LOWB.GE. BUBS IGO  TO  6200 

00003850 

C 

CHECK  IF  CURRENT  SOLUTION  SATISFIES  CAPACITY 

00003880 

C 

CONSTRAINTS 

00003890 

0202 

4200 

IF  (IUNCAP.EC.il  GO  TO  4420 

00003900 

0283 

4210 

00  4400  K=1,P 

00003910 

0284 

NSUMDsO 

00003920 

0205 

00  4300  J* 1 »N 

00003930 

C2Ue 

IX  =SCLX ( J 1 

00003950 

0287 

NSUMD  =NSUMD*D ( 1 X , J ,  X) 

00003960 

02  08 

4300 

CONTINUE 

00003970 

0269 

IF  (  ISTEP.Nfc.2l  GO  TO  4320 

00003980 

0290 

PRINT  4310*  K, NSUMD 

00003990 
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FORMAT  ('O',  'K  »  NSUMO 1 «  2 1 10) 

IF(N$UMO.LE.S(K)t  GO  TO  4400 
GO  TO  5100 
CONTINUE 

********** ****CCMPUT£  UPPER  BOUND  UPB  IF  CAPACITY  CONSTRAINTS 
ARE  SATISFIED. 

UPB*SUM  OF  A(I,J)*F1XED  COST  FCU3  BASED  ON 
SOLUTION  VECTOR  SOLX(J) 

VECTOR  OF  FACILITIES  FOR  UPPER  BOUND  FUB(K>  MAS 
VALUES  1  OR  0  BASED  ON  FACITILY  USED  OR  OTHERWISE 

DO  4450  K=1,P 
FUB(K 1*0 
CONTINUE 
NSUMA=0 
FCUB*0 

DO  4650  J*1,N 
IX=SOLX (J) 

NSUMA*NSUM A4A ( IX  t  J I 
00  4600  K*1,P 

IF  (E(  IX,K) .EQ.O)  GO  TO  4600 
IF ( FUB (KI.EQ.l)  GO  TO  4600 
FUB(KI*1 
FCU8*FCU8+6 (K ) 

CONTINUE 

CONTINUE 

IF  (ISTEP.NE.2I  GO  TO  4700 
PRINT  4660,  (FUBIK) ,K=1,P) 

FORMA  T('0'«'( FUB ( K  J ,K=  1  ,  P  )  ',  2014/16X, 20141 

UPB*NSUMA*FCU8 

IF  (1STEP.E0.0I  GO  TO  4750 

PRINT  4710,  NSUMA,  FCUB ,  UPB,  BUB,  BUBS 

FORMAT  ('O',  'NSUM A ,  FCUB,  UPB,  BUB,  BUBS  *,2110,  .'F15.4J 
COMPARE  UPPER  BOUND  WITH  BEST  UPPER  BOUND 
IF  UPB  IS  LESS  THAN  BUB,  SET  IT  AS  BUB  ANO 
NOTE  THE  SOLUTION 
IF  (UPB. GE. BUB)  GO  TO  5100 
BUB*UPB 

BUBS*  BUB/  (1.0*EPS) 

IBNOO*NOD 

PRINT  4780,  IBNOO,  BUB,  BUBS 

FORMAT  ('O',  • IBNQD,  BUB,  BUBS',  110,  2F15.4) 

□0  4800  J*1,N 
BSOLX ( J ) *SOLX ( J ) 

00  4850  K*1 ,P 
BSOLY (K)»FU8(K) 

•  •************COMP A  RE  LOWB  WITH  BUBS.  IF  LOWB  IS  GREATER 
THAN  OR  EQUAL  TO  BUBS,  THEN  8ACXTP. ACX 
IF  (LOWB .GE . BU8S )G0  TO  6200 

FIX(J)  VALUES  ARE  1  FOR  EACH  J,  THEN  BACKTRACK 
IF  (LQ1.E0.N)  GO  TO  6200 

*••*••******•* APPLY  THE  BOUNCING  RULE************«************ 

IF  THE  DIFFERENCE  BETWEEN  C(I,J>  AND  MINC(J)  IS 
GREATER  THAN  THE  DIFFERENCE  BETUtEN  BUBS  AND 
LOWB,  THEN  C* ( I  ,  J  )  *2 

♦********R****APPLY  BRANCHING  RULE  AND  FIND  BRl,  THE  NEXT 
VARIABLE  FDR  LEFT  BRANCHING. 

FINO  NMINCUI,  THE  NEXT  HIGHER  VALUE  THAN  MINC(J) 
ANO  OIFBR(J),  THE  DIFFERENCE  BETWEEN  THEM. 


00004000 
00004010 
00004020 
00004030 
00004040 
00004050 
00004060 
00004070 
00004080 
00004090 
00004100 
00004110 
00004120 
00004130 
00004140 
OOC  4150 
00004170 
000041B0 
00004190 
00004200 
00004210 
00004220 
00004230 
00004240 
00004250 
00004260 
00004270 
00004280 
00004290 
00004300 
00004310 
00004320 
00004330 
00004340 
00004350 
00004360 
00004370 
00004380 
00004385 
00004386 
00004388 
00004390 
00004400 
00004410 
00004420 
00004430 
00004440 
00004450 
00004480 
00004500 
00004510 
00004515 
00004520 
00004525 
00004530 
00004540 
00004550 
00004555 
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0329 

0330 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

033S 

0339 

0340 

0341 

0342 

0343 

0344 

0345 

0346 

0347 

0343 

0349 

0350 

0351 

0352 

0353 

0354 

0355 

0356 

0357 

0358 

0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

0368 

0369 
0370 
0371 
0372 
0373 
03  74 
0375 
0376 
0377 
0378 
03  79 
0330 
0381 
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5300 

5320 

5330 


D80UND*BU8S“L0W8 
5200  00  5250  J*l,N 

NH1NC I J 1=0.0 
01 F8R I J  >  =0.0 
5250  CONTINUE 

00  5600  J«1,N 

C  SKIP  TO  NEXT  J  IF  FIX(J»=l 

IF  (FIX(J).EQ.l)  GO  TO  5600 
«-K*0 

c  "l  SKIP  Cd.JJ  IF  CXd,J»*2  C  HOVE  TO  NEXT  I 

IF  ( CXt I« J) .EQ.2 )  GO  TO  5350 
IF  d.EC.SOLX(J) »  GO  TO  5350 
IF  ((CUfJ)-HINCUn  .GT.D90UND)  GO  TO  5330 
IF  HI-UO.EQ.l)  GO  TO  5320 
IF  (C ( 1 , J ) .GE.NMINC ( J I )  GO  TO  5400 
NM INC (4)*C4I#JI 
GO  TO  5400 
CX«I ,J»=2 
NSX-NSX+1 

STX(NSXI=-d*100*J 1-1000000 
KT2(Jl«KT24Jl-»l 

IF(KT2U».LT.<M-II»  GO  TO  5350 
INDI=SOLX( J) 

CX(INDI,J»=1 
N'.X=N$X*l 

STX ( NSX )  =  4 INOI* 100*J  >  *1000000 
FIX(JJ=l 
LQ1=LQI+1 
FIXI ( J )=1N0I 
GO  TO  SoOO 
LK=LK*l 
l-l*l 

IFII.LE.H)  GO  TO  5300 
DIF8RU>=NMINC(J>-MINC(J> 

CONTINUE 

IF  IISTEP.NE.2I  GO  TO  5650 
DO  5620  1*1, M 

PRINT  290,  I,tCXd,J»,J*l.NI 
CONTINUE 

PRINT  297,  (FIX(J)  * J* 1 «N I 

IF  FIX(JI=1  FOR  ALL  J,  THEN  BACKTRACK, 

IF  (LQl.EQ.Nl  GO  TO  6200  .....  .. 

FIND  HAXDI F  »  THE  MAXIMUM  DIF-ERENCE  DIFBR(J) 

LF=0 

DO  5800  J=1 ,N 

IF  (FIXIJt.EO.ll  GO  TO  5690 
IF  4  I J-LF ) .EQ.l )  GO  TO  5660 
IF  (DIFBRUI.lt. MAXDIF)  GO  TO  5800 
5660  MAXOI F=DI F8R ( J I 


5350 

5400 

5500 

5600 


5620 


5650 


LJ=J 

GO  TO  5800 
5690  LF*LF*l 
5800  CONTINUE 

IF  !ISTEP.Nt.2)  GO  TO  5840 
DO  5820  J*1,N 

IF  (FlXUl  .EQ.l  »  GO  TO  5820 

PRINT  5810,  J,  NMINCUJ,  MINCIJI,  OIFBR(J) 


00004568 

00004570 

00004580 

00004590 

00004600 

00004610 

00004620 

00004630 

00004640 

00004650 

00004670 

00004680 

00004690 

00004700 

00004710' 

00004720 

00004730 

00004735 

00004740 

00004742 

00004745 

00004747 

00004750 

00004752 

00004755 

00004758 

00004760 

00004762 

00004764 

00004766 

00004768 

00004770 

00004775 

00004780 

00004785 

00004790 

00004795 

00004820 

00004830 

00004850 

00004860 

00004880 

00004890 

00004900 

00004905 

00004910 

00004915 

00004920 

00004925 

00004930 

00004935 

00004940 

00004943 

00004946 

00004950 

00004953 

00004956 

00004960 
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0383 

5810 

FORMAT  ( *0*, *J,NMINC(JI,MINC( JI.DIFBRIJ) I5.3F15.4I 

00004963 

0384 

5820 

CONTINUE 

00004966 

C 

00004970 

0385 

5840 

DO  5900  J»1,N 

00004980 

03b6 

IF  (J.NE.LJ)  GO  TO  5900 

00004990 

0387 

5850 

BR1=S0LX( JI*100*J 

00005000 

0388 

IF  (ISTEP. EQ. 01  GO  TO  6020 

00005010 

0389 

PRINT  5880,  BR1 

00005020 

0390 

5880 

FORMAT ( *0 ' » •  8R l • , I 10 ) 

00005030 

0391 

GO  TO  6020 

00005040 

0392 

5900 

CONTINUE 

00005050 

C 

00005060 

C 

NSX  REPRESENTS  THE  NUMBER  OF  VARIABLES  IN  STXIINSI 

00005070 

0393 

6020 

NSX=NSX*1 

00005090 

0394 

6040 

STX 1 NSX ) *BR1 

00005100 

0395 

IF  IISTEP.NE.2)  GO  TO  6100 

00005150 

0396 

PRINT  6088,  ( STX ( INS ) ,  INS=1,NSX! 

00005160 

0397 

o088 

FORMAT  1 • 0* , •  STXIINSI1,  10110,  1221/,  10X,10I10II 

00005170 

C 

00005220 

C 

CAPACITY  RULE 

00005230 

0398 

6100 

NOD=NOD+ l 

00005240 

0399 

6110 

IF  (ET.EQ.O.OI  GO  TO  6150 

00005242 

0400 

IFIlNSET.EO.il  GO  TO  6147 

00005244 

0401 

IF  (1NET.E0.1I  GO  TO  6150 

00005246 

0402 

CALL  TIMETTINTI 

00005248 

0403 

EL TN= ( INT-I TO  I *26.046-6 

00005250 

0404 

IF  IELTN.LT.ETI  GO  TO  6150 

00005253 

0405 

6120 

PRINT  6125,  NOO,  ELTN,  BUB,  BUBS,  IBNOD 

00005256 

0406 

6125 

FORMAT  ( *0  • ,  ‘WAS  AT  NOOE',16,  •  AT  ELAPSED  TIME  =*,  F10.4, 

00005260 

1 

•  SECQNOS.  '  ,/lX,  •  BUBs'.FlS.*,  •  BUB S=  •  , F 1 5 . 4 , 

00005263 

2 

•  AT  N0DE=«,17I 

00005266 

0407 

I BUBX8UB 

00005267 

0408 

IF  (IBUB.EQ. 99999991  GO  TO  6146 

00005268 

0409 

6130 

PRINT  6135,  IBSOLXI J),J=1,N1 

00005270 

0410 

0135 

FORMAT ("O',  •SOLUTION  CORRESPONDING  TO  BUB  IS',  //IX, 

00005273 

1 

•IBSOLXIJI,  J  =  l,N I  •, 1018 ,3 I/18X, 1018  I  1 

00005276 

0411 

6140 

PRINT  6145,  (BSOLYIKI,  K=l,Pl 

00005280 

0412 

6145 

F0RMAT(/1X,MBS0LY(KI,  K  =  1,P|  *,10I8,2(/iaX,  1018  1 1 

00005290 

0413 

6146 

I NE  T* 1 

00005292 

0414 

IN IS=  ISTEP 

00005294 

0415 

I NSET=1 

00005296 

0416 

ISTEP=2 

00005298 

0417 

GO  TO  6150 

00005300 

0418 

6147 

ISTEP=INIS 

00005302 

0419 

1NSET=0 

00005304 

0420 

6150 

GO  TO  272 

00005306 

C 

00005308 

0421 

6200 

IF  (NSX.EQ.OI  GO  TO  8100 

00005310 

0422 

6220 

IF  (  IABSI STXINSXI I.GT. 10000001  GO  TO  6500 

00005320 

0423 

62  50 

BR0=STX ( NSX 1 

00005330 

0424 

62  70 

S  TX (NSX M-BRO-tOCOOOO 

00005340 

0425 

IF  (ISTEP. EO. 01  GO  TO  6308 

00005390 

0426 

PRINT  6305,  BRO 

00005400 

0427 

6305 

FOR  MA  T ( *0 • ,  *BR0  *,110) 

00005410 

0428 

6308 

IF  llSTfP.NE.2l  GO  TO  6330 

00005*20 

0429 

PRINT  6088,  (STXIINSI,  INS=I,NSXI 

00005430 

C 

00005490 

C 

CAPACITY  rule 

00003500 

94 
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0430 

6330 

NODsNDOU 

00005510 

0431 

6410 

IF  (ET. £0.0.0)  GO  TO  6450 

00005512 

l 

04  32 

IF  UNSET.  EQ.l)  GO  TO  6445 

00005516 

0433 

IF  (INET.EC.l)  GO  TQ  6450 

00005518 

0434 

CALL  TtHETClNTI 

00005520 

t. 

04  35 

ELTN«(INT-lT0)*26.04E-6 

00005523 

0436 

IF  (ELTN.LT.tT)  GO  TO  6450 

00005526 

0437 

6420 

PRINT  6125,  NOD,  ELTN,  BUB,  BUBS,  IBNOO 

00005528 

0 

0438 

IBUB*BUB 

00005530 

0439 

IF  UBUB.EO. 9999999)  GO  TO  6442 

00005532 

0440 

6430 

PRINT  6135,  (BSOLX(J) , J-l ,N) 

00005533 

© 

0441 

6440 

PRINT  6145,  (BSOLY(K),  K=1,P) 

00005536 

0442 

6442 

INE  T= 1 

00005538 

0443 

INIS*IST£P 

00005540 

0 

0444 

INSET  =  1 

00005542 

0445 

TSTEP*2 

00005544 

0446 

GO  TO  6450 

00005546 

o 

0447 

6445 

ISTEP=IN1S 

00005548 

0448 

INSET  =0 

00005550 

0449 

6450 

GO  TO  210 

00005552 

c 

0450 

6500 

IF  C  STX(NSX) .GT. 1000000)  GO  TO  6520 

00005555 

0451 

LX=-STX(NSX 1-1000000 

00005560 

0452 

I X  =LX/100 

00005570 

0453 

JX*LX-IX*100 

00005580 

0454 

CXtIX,JX>=0 

00005590 

0455 

KT2(JXI=KT2<JXI-1 

00005595 

r 

04  56 

GO  TO  6550 

00005600 

0457 

6520 

LX*  STX(NSX)— 1000000 

00005610 

0**58 

IX=LX/100 

00005620 

■ .  . 

0459 

JX=LX-IX*100 

00005630 

0460 

CXIIX,JX)*0 

00005640 

0461 

F.X( JX)=0 

00005650 

sj 

0462 

LQ1=L01-1 

00005660 

0463 

6550 

NSX=NSX-1 

00005690 

0464 

GO  TO  6200 

00005700 

C 

00005730 

0465 

8100 

1BUB=BUB 

00005740 

0466 

CALL  TIMET ( I Tl ) 

00005750 

0 

0467 

ELT1=< ITl-IT0)*26.04E-6 

00005760 

0468 

PRINT  8105,  ELT1 

00005770 

0469 

8105 

FORMAT  C'0',///1X,  'ELAPSED  TIME  IN  SECONDS-',  F15. 

8) 

00005780 

0 

0470 

PRINT  8120,  NOO 

00005790 

0471 

8120 

FORMAT  (  'O', 'TOTAL  NUMBER  OF  NODES  EXPLORED  *»,I3) 

00005800 

04  72 

IF  ( IBUB.EQ. 9999999)  GO  TO  8350 

00005810 

c 

0473 

8130 

PRINT  8150 

00005820 

0474 

8150 

FORMAT  I 'O',  'NOTE:  1.  FOLLOWING  XCI,J) 

VARIABLES 

SHOW  DESIGN', 

00005830 

l  •  I  TO  WHICH  ACTIVITY  J  IS  ASSIGNED  FOR  J*1 

TO  N. ', 

00005840 

© 

2  /7X »  '2.  IF  EPSILON  EPS  WAS  ASSIGNED  A  POSITIVE', 

00005850 

3  '  (NON-ZERO)  VALUE,  THE  SOLUTION 

MAY  BE  SUBOPTIMAL.',/! 
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